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FOREWORD 


This  document  is  the  final  report  for  U.S.  Army  Research  Office  contract 
number  DAAG29-85-K-0054,  which  corresponds  to  the  proposal  entitled  Rapid 
Thermal  Processing  of  III-V  Compound  Semiconductors  With  Appli¬ 
cation  to  the  Fabrication  of  Microwave  Devices.  This  research  is  focused 
on  establishing  better  techniques  for  surface  doping  of  III-V  semiconductors.  Por¬ 
tions  of  this  work  were  carried  out  in  conjunction  with  our  research  on  Limited 
Reaction  Processing  and  MOCVD  of  compound  semiconductors,  under  a  previous 
ARO  contract,  DAAG29-85-K-0237. 

The  GaAs  doping  techniques  investigated  in  this  work  include  Si  surface 
diffusions  using  a  solid  source,  ion  implantation  followed  by  Rapid  Thermal  An¬ 
nealing,  and  both  solid  source  and  vapor  phase  Zn  diffusions.  New  models  for  Si 
diffusion/electrical  activation,  and  Zn  diffusion  have  been  developed  and  imple¬ 
mented  in  computer  simulations.  Using  these  rapid  thermal  diffusion  techniques, 
we  have  demonstrated  that  shallow  n+  and  p+  layers  can  be  produced  reliably  while 
retaining  excellent  surface  morphology.  In  addition,  an  improved  understanding 
of  the  annealing  of  defects  and  electrical  activation  of  dopants  in  ion  implanted, 
amorphized  GaAs  has  been  achieved. 

This  report  constitutes  a  brief  summary  of  research  performed  under  the 
above  contract.  It  begins  in  Part  1  with  a  list  of  publications  associated  with  this 
research.  A  short  synopsis  of  7  publications  which  describe  our  findings  is  followed 
by  Part  3,  which  contains  reprints  of  these  publications.  The  final  section  (Part  4) 
contains  a  list  of  participating  scientific  personnel. 
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Part  1 

List  of  Publications 


The  following  is  a  list  of  7  publications  which  are  summarized  in  Part  2  of 
this  report.  Reprints  of  these  articles  appear  in  Part  3  in  the  following  order: 

1.  Mark  E.  Greiner  and  James  F.  Gibbons,  “Diffusion  and  electrical  prop¬ 
erties  of  silicon-doped  gallium  arsenide'’,  J.  Appl.  Phys.  57,  15  June 
1985,  pp.  5181-5187. 

2.  W.G.  Opyd  and  J.F.  Gibbons,  “Regrowth  of  amorphized  compound 
semiconductors",  Mat.  Re3.  Soc.  Symp.  Proc.  Vol.  45,  (Materials 
Research  Society,  Pittsburg,  1985),  pp.  273-277. 

3.  W.G.  Opyd,  J.F.  Gibbons,  J.C.  Bravman,  and  M.A.  Parker,  “Damage 
calculation  and  measurement  for  GaAs  amorphized  by  Si  implantation" . 

Appl.  Phys.  Lett.  49,  13  October  1986,  pp.  974-976. 

4.  W.G.  Opyd,  J.F.  Gibbons,  and  A.J.  Mardinly,  “Precipitation  of  Impu¬ 
rities  in  GaAs  Amorphized  by  Ion  Implantation",  submitted  to  Appl. 

Phys.  Lett.,  May  1988. 

5.  S.  Reynolds,  D.W.  Vook,  and  J.F,  Gibbons,  “Limited  reaction  process¬ 
ing:  Growth  of  III-V  epitaxial  layers  by  rapid  thermal  metalogranic 
chemical  vapor  deposition”,  Appl.  Phys.  Lett.  49,  22  Dec.  1986,  pp. 
1720-1722. 

6.  S.  Reynolds,  D.W.  Vook,  W.G.  Opyd,  and  J.F.  Gibbons,  “Rapid  ther¬ 
mal  annealing  of  Si-implanted  GaAs  with  trimethylarsenic  overpres¬ 
sure”,  Appl.  Phys.  Lett.  51,  21  Dec.  1987,  pp.  916-918. 

7.  S.  Reynolds,  D.W.  Vook,  and  J.F.  Gibbons, “Open-tube  Zn  diffusion  in 
GaAs  using  diethylzinc  and  trimethylarsenic:  Experiment  and  model", 

J.  Appl.  Phys.  63,  15  Feb.  1988,  pp.  1052-1059. 
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Summary  of  the  Research 


Rapid  Thermal  Processing  (RTP)  offers  several  important  advantages  over 
furnace  annealing  of  GaAs.  In  fact,  because  of  difficulties  associated  with  surface 
decomposition,  the  arguments  in  favor  of  employing  RTP  are  often  more  compelling 
for  compound  semiconductors  than  for  silicon.  For  a  bare  surface,  arsenic  and  gal¬ 
lium  loss  rates  are  very  high  (  1016  /cm2  —  sec.  at  800  °C)  at  temperatures  required 
to  anneal  implant  damage  or  perform  dopant  diffusion  in  GaAs.  Hence,  limiting 
the  time  at  temperature  is  important  for  GaAs  processing.  In  addition,  because  of 
the  system  design,  it  is  usually  easier  to  control  the  ambient  during  RTP  compared 
to  furnace  processing.  Using  RTP,  we  have  also  been  able  to  study  fundamental 
processes  such  as  diffusion  and  implant  damage  annealing  in  a  time/temperature 
regime  not  accessible  with  conventional  techniques. 

2.1  Si  Diffusion  in  GaAs 

Applications  of  thin,  heavily  Si-doped  layers  include  source  and  drain  con¬ 
tacts  for  GaAs  MESFETS  .  Work  performed  as  part  of  a  previous  ARO  contract 
included  a  study  of  RTP  employed  to  diffuse  Si  into  GaAs  from  a  thin  elemental 
source.  Si3N4  and  SiC>2  encapsulants  were  used  in  these  experiments.  Diffusions 
were  performed  at  temperatures  as  high  as  1050  CC  for  times  as  short  as  3  sec. 
Si  was  observed  to  diffuse  into  GaAs  at  depths  up  to  0.3  /im  for  SiC>2  capped 
samples.  Under  otherwise  identical  conditions,  no  diffusion  was  detected  for  Si3N4 
encapsulated  GaAs. 

In  the  early  part  of  the  present  contract  period,  we  developed  a  complete 
diffusion  model  for  Si  in  GaAs  to  explain  the  measured  RTA  data.  This  model  is 
described  in  Paper  #  1  in  Section  3.  The  Si  diffusion  mechanism  is  based  upon 
the  formation  of  nearest-neighbor  donor- acceptor  pairs.  The  diffusion  front  is  very 
steep  due  to  pair  dissociation  at  low  dopant  concentrations.  Several  experiments 


3 


\ 


* 


which  support  the  application  of  this  mechanism  to  high  concentration  Si  diffusion 
in  GaAs  are  discussed.  A  compensation  mechanism  for  amphoteric  dopants  is  de¬ 
veloped  as  well.  These  models  are  potentially  applicable  to  all  amphoteric  dopants 
in  III-V  compounds  where  dopant  pairing  is  likely  to  dominate  the  diffusion  pro¬ 
cess. 


2.2  Solid  phase  epitaxy  of  compound  semiconductors 

Paper  #  2  describes  the  initial  motivation0  for  an  investigation  of  solid  phase 
epitaxy  (SPE)  in  compound  semiconductors,  which  include  the  possibility  of  lower 
temperature  anneals,  as  in  the  case  of  silicon,  which  would  result  in  reduced  ma¬ 
terial  decomposition  and  impurity  diffusion.  The  thresholds  for  amorphization 
of  InSb  and  GaAs  are  reported  based  upon  correlations  of  Boltzmann  transport 
calculations  with  experimental  observations. 

Paper  #  3  presents  a  study  of  the  stages  of  defect  annealing  following  the 
SPE  of  an  amorphous  layer.  Material  is  analyzed  over  a  range  of  annealing  temper¬ 
atures  by  transmission  electron  microscopy  (TEM)  and  Rutherford  backscattering 
spectrometry.  Though  SPE  occurs  at  low  temperatures,  residual  defect  annealing 
is  noted  to  require  higher  temperatures. 

Paper  #  4  describes  an  investigation  of  the  problem  of  failure  to  achieve 
electrical  activation  of  implanted  impurities  even  after  high  temperature  defect  an¬ 
neals.  Photoluminescence,  TEM,  and  x-ray  microanalysis  are  used  to  show  that 
the  lack  of  electrical  activation  is  associated  with  the  precipitation  of  stable,  elec¬ 
trically  neutral  impurity  complexes.  This  result  is  shown  to  be  similar  to  that 
observed  for  low  temperature  MBE  and  demonstrates  a  technique  for  precipitating 
impurities  in  GaAs. 

As  part  of  a  companion  contract,  we  have  developed  a  III-V,  Rapid  Thermal 
MOCVD  system,  as  described  in  Paper  #  5.  Besides  epitaxial  growth,  we  were  also 
able  to  use  this  system  to  study  Rapid  Thermal  Annealing  of  Si  implants  in  GaAs 
in  a  trimethylarsenic  overpressure.  As  described  in  Paper  #  6,  we  have  compared 
Si  implant  activation  efficienty  and  surface  degredation  for  arsenic  ambient  and 
proximity  capped  anneals.  The  arsenic  ambient  gives  consistenly  higher  implant 
activation  with  better  surface  morphology. 
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2.3  Zn  Diffusion  in  GaAs 


Zn  diffusion  in  GaAs  is  an  important  processing  step  in  the  fabrication  of 
diodes,  junction  field-effect  transistors,  and  bipolar  transistors.  Using  conven¬ 
tional  sealed  ampoule  methods,  control  over  the  diffusion  process  is  often  difficult, 
and  poor  surface  morphology  and  variable  junction  depths  can  limit  device  perfor¬ 
mance.  Using  the  organometallic  sources  diethylzinc  (DEZn)  and  trimethylarsenic 
(TMAs),  we  have  characterized  the  diffusion  of  Zn  into  GaAs,  as  described  in  Pa¬ 
per  #  7.  This  method  produces  surface  hole  concentrations  in  excess  of  1020  cm-3, 
with  good  control  of  junction  depths  as  shallow  as  0.1  pm.  A  key  advantage  of 
this  method  is  the  smooth  surface  morphologies  obtained. 

As  discussed  in  Paper  #  7,  diffusions  were  carried  out  in  a  hot-wall  quartz 
reaction  chamber  at  atmospheric  pressure.  The  gas  delivery  system  is  essentially 
the  same  one  which  is  described  in  Papers  5  and  6.  Using  our  MOCVD  reactor 
to  perform  vapor  phase  diffusions  allows  us  to  accurately  meter  the  quantities  of 
reactants  delivered  to  the  wafer.  We  were  able  to  independently  control  the  time, 
temperature,  and  partial  pressures  of  As  and  Zn,  as  well  as  the  total  pressure, 
which  can  not  be  done  using  the  conventional  ampoule  method.  In  total,  over  40 
diffusion  samples  were  analyzed  by  SIMS  and  Hall  effect  measuiements. 

Paper  #  7  describes  the  results  in  detail.  These  diffusions  have  a  com¬ 
plex  double  profile  which  is  not  predicted  by  the  widely  accepted  interstitial- 
substitutional  model  for  Zn  diffusion.  We  suggest  that  this  double  profile  may 
be  caused  by  the  existence  of  two  mobile  Zn  species  which  dominate  in  different 
regions  according  to  the  position  of  the  local  Fermi  level.  We  have  implemented  a 
computer  simulation  which  gives  an  approximate  solution  based  upon  this  model, 
and  provides  a  good  fit  to  both  our  data,  and  that  of  several  other  authors.  Our 
model  shows  that  a  double  diffusion  profile  does  not  necessarily  imply  vacancy  (or 
interstitial)  nonequilibrium. 
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Reprints  of  Publications 


This  sections  contains  reprints  of  the  articles  listed  in  Part  1  and  described 
in  Part  2. 
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Payer  #1 

Diffusion  and  electrical  properties  of  silicon-doped  gallium  arsenide 

Mark  E.  Greiner"  and  James  F.  Gibbons 

Stanford  Electronics  Laboratory ,  Stanford  University,  Stanford,  California  94305 
(Received  9  October  1984;  accepted  for  publication  4  January  1985) 

The  amphoteric  nature  of  silicon  in  gallium  arsenide  is  used  to  develop  diffusion  and  electrical 
compensation  mechanisms.  The  diffusion  mechanism  is  based  on  the  formation  and  diffusion  of 
nearest-neighbor  donor-acceptor  pairs.  General  solutions  are  presented  that  predict  abrupt 
diffusion  fronts  for  a  wide  range  of  pairing  conditions.  Experiments  support  the  application  of  this 
mechanism  to  Si  diffusion  in  GaAs  at  high  concentrations.  A  compensation  mechanism  for 
amphoteric  dopants  is  developed  as  well.  The  compensation  process  is  driven  primarily  by  the 
free-electron  concentration.  Nearly  complete  compensation  is  predicted  for  large  dopant 
concentrations. 


I.  INTRODUCTION 

Except  for  a  few  attempts,1'1  Si  diffusion  has  not  been 
used  to  form  n-type  layers  in  GaAs.  Because  of  this,  a  well- 
developed  technology  for  Si  diffusion  does  not  exist  at  this 
time.  The  current  processes  for  forming  n-t ype  layers  in 
GaAs  involve  the  growth  of  epitaxial  layers  and  the  implan¬ 
tation  of  n-type  dopants  into  semi-insulating  GaAs.  In  the 
latte  case,  studies  were  done4-7  on  the  Si  diffusion  which 
occurs  during  subsequent  annealing  of  these  ion-implanted 
layers.  The  concentrations  of  implanted  Si  in  these  studies 
are  limited  to  those  used  in  GaAs  field  effect  transistor  chan¬ 
nels  and  are  generally  in  the  10l7/cm3  range.  This  study  fo¬ 
cuses  on  Si  in  GaAs  diffused  from  the  surface  under  condi¬ 
tions  that  result  in  heavily  doped  n-type  layers.  A  diffusion 
model  and  a  compensation  model  are  developed  for  ampho¬ 
teric  dopants  and  applied  to  Si  doping  in  GaAs. 

Before  presenting  the  results  from  this  work,  two  sets  of 
published  diffusion  data  will  be  reviewed.  The  first  set  con¬ 
tains  the  early  diffusion  work  by  Vieland1  and  Antell,2  who 
use  Si  vapor  and  solid  sources,  respectively,  to  form  diffused 
layers,  and  recent  work  by  the  present  authors.3  The  second 
set  of  data  to  be  reviewed  deals  with  diffusion  of  low-dose, 
ion-implanted  Si.  These  results  will  be  useful  as  a  compari¬ 
son  to  the  diffusion  effects  observed  at  high  concentrations. 

II.  REVIEW  OF  Si  DIFFUSION 
A.  High  SI  concentrations 

Both  Vieland  and  Antell  reported  the  effect  of  arsenic 
pressures  on  diffusion  rates  in  GaAs.  Vieland  diffused  Si  into 
GaAs  at  1 1 50  *C  and  relied  on  the  vapor  pressure  of  Si  over  a 
solid  Si  source  ( ~  10~6  Torr)  to  provide  the  dopant.  In  An- 
tell’s  work.  Si  was  diffused  from  a  3000-A  Si  film  sputtered 
directly  onto  a  p- type  GaAs  sample.  All  anneals  took  place 
in  evacuated  ampoules  containing  solid  arsenic  to  produce 
the  desired  arsenic  overpressures.  For  all  three  temperatures 
studied,  the  junction  depth  increased  rapidly  with  increasing 
arsenic  overpressure  until  a  knee  was  reached  at  some  pres¬ 
sure  PK .  Above  this  pressure,  the  junction  depth  became 
independent  of  pressure.  Without  a  detailed  knowledge  of 


the  GaAsSi  ternary  diagram,  it  is  difficult  to  interpret  this 
knee. 

The  arsenic  overpressure  is  predicted  to  vary  the  con¬ 
centration  of  gallium  and  arsenic  vacancies  by  the  interac¬ 
tion  between  the  arsenic  vapor  and  the  GaAs  surface.  The 
magnitude  of  these  variations  is  expressed  in  the  following 
equations1: 

[F0.  ]=*,./»£  (1) 

and 

[r*]=k2.pz!'\  (2) 

By  comparing  this  vacancy  dependence  on  As  overpressure 
with  the  experiments  of  Vieland  and  Antell  above,  we  see 
that  the  diffusion  rate  of  Si  increases  with  the  surface  con¬ 
centration  of  gallium  vacancies. 

In  addition  to  the  junction  depth,  Antell  measured  a 
free-carrier  concentration  of  5  X  10'Vcm3  for  a  1000  °C  dif¬ 
fusion.  This  concentration  was  constant  throughout  the  dif¬ 
fused  layer.  The  surface  concentration  was  estimated  to  be 
greater  than  10‘°/cm3  showing  a  high  level  of  electrical  com¬ 
pensation  in  these  diffused  layers. 

In  a  recent  study  by  the  present  authors3  rapid  thermal 
processing  was  used  to  diffuse  Si  into  GaAs  from  a  thin  ele¬ 
mental  source.  The  sample  structure  is  shown  in  Fig.  1. 
SijN«  and  Si02  encapsulants  were  used  here  rather  than  As 
overpressures.  Diffusions  were  performed  at  temperatures 
as  high  as  1 050  °C  for  times  as  short  as  3  sec.  Si  diffusion  into 
GaAs  up  to  0. 3  p  m  was  observed  for  Si02  encapsulated  sam¬ 
ples.  No  diffusion  was  detected  for  SijN4  encapsulated  sam¬ 
ples  for  the  short  time  used  in  this  process. 


Si02  Of  S13N4 
100  A  Si 


GaAs 


*'  Current  address.  Cincinnati  Electronics  Corporation,  2630Glendale-Mil-  FIG  1  Sample  structure  used  to  diffuse  Si  into  GaAs  during  rapid  thermal 
ford  Road,  Cincinnati,  Ohio  45241.  processing 
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Secondary  ion  mass  spectroscopy  (SIMS)  was  used  to 
measure  the  St  diffusion  profiles.  The  resulting  profiles  indi¬ 
cated  a  concentration-dependent  diffusion  coefficient.  The 
Si  diffusion  front  was  very  steep,  similar  to  the  diffusion  pro¬ 
files  produced  by  zinc  diffusion  in  GaAs.  Although  both  Zn 
and  Si  exhibit  these  abrupt  diffusion  fronts,  they  behave  dif¬ 
ferently  under  As  overpressures.  Zn  is  believed  to  diffuse  by 
an  interstitial-substitutional  mechanism. K  Increasing  the  As 
overpressure  retards  the  Zn  diffusion  process,9  in  agreement 
with  the  interstitial-substitutional  mechanism.  In  contrast, 
Si  diffusion  is  enhanced  under  As  overpressures. 12  Also,  the 
activation  energy  for  the  Si  diffusion  process  at  high  concen¬ 
trations  was  measured  to  be  2.S  eV,3  which  is  in  the  range 
associated  with  substitutional  diffusion  in  GaAs. 

A  substitutional  mechanism  was  proposed3  for  Si  diffu¬ 
sion  in  GaAs  based  upon  the  formation  and  diffusion  of  Si 
nearest-neighbor  pairs.  This  model  assumes  that  Si  pairs 
dominate  the  diffusion  flux  at  high  Si  concentrations.  Since 
the  pair  concentration  is  a  function  of  the  total  Si  concentra¬ 
tion,  a  concentration-dependent  diffusion  coefficient  is  pre¬ 
dicted  by  this  mechanism.  This  model  provides  a  good  fit  to 
the  profiles  produced  by  Si  diffusion  during  rapid  thermal 
processing. 

8.  Low  Si  concentrations 

There  have  been  many  reports  on  the  electrical  profiles 
produced  by  implanting  Si+  into  semi-insulating  GaAs.4'7 
After  the  implant,  the  silicon  is  made  electrically  active  by  a 
postimplant  anneal.  These  anneals  use  some  method  de¬ 
signed  to  eliminate  surface  erosion  due  to  preferential  As 
evaporation.  The  resulting  electrical  profiles  depend  on  the 
method  used  during  the  postimplantation  anneal.  As  a  re¬ 
sult,  reported  diffusion  coefficients  vary  over  a  wide  range. 
The  results  from  several  authors  are  summarized  in  Table  I. 

No  two  sets  of  experiments  were  found  with  identical 
implant  and  annealing  conditions.  The  results  in  Table  I  cov¬ 
er  the  range  of  techniques  commonly  used  in  furnace  anneal¬ 
ing.  In  all  cases  when  diffusion  was  measurable,  the  as-im¬ 
planted  profiles  were  simply  broadened  while  maintaining 
their  as-implanted  shape.  The  retention  of  the  as-implanted 
shape  under  all  annealing  conditions  implies  there  is  a  com¬ 
mon  diffusion  mechanism  in  which  the  diffusion  coefficient 
of  Si  is  not  a  function  of  concentration  over  the  range  of 
concentrations  studied.  These  concentrations  are  less  than 
lO'VcmL  However,  concentration-dependent  diffusion  was 
observed  at  Si  concentrations  between  10'Vcm  and  10:o/ 


The  value  of  the  diffusion  coefficient  for  these  low  Si 
concentrations  depends  upon  the  specific  annealing  condi¬ 
tions.  In  a  study  by  Kasahara  and  Watanabe."’  the  diffusion 
of  low-dose  Si’’’  implants  in  unencapsulated  GaAs  was  mea¬ 
sured  after  annealing  at  900 “C  under  various  partial  pres¬ 
sures  of  AsH,.  The  diffusion  rate  was  found  to  increase  with 
the  fourth  root  of  the  AsH,  partial  pressure  for  pressures 
greater  than  0.2  Torr.  Two  of  these  values  are  listed  in  Table 
I.  This  dependence  is  similar  to  that  obtained  by  Vieland  and 
Antell  for  high-concentration  diffusion  under  As  overpres¬ 
sures. 

The  diffusion  data  for  Onuma  et  a l.1  in  Table  I  for  en¬ 
capsulated  anneals  show  a  dependence  on  the  type  of  encap- 
sulant.  SiO:  encapsulated  anneals  produced  Si  diffusion 
rates  20  times  greater  than  anneals  using  SijN4  encapsula¬ 
tion.  It  was  suggested7  that  this  enhanced  diffusion  under 
SiO,  encapsulation  was  related  to  the  excessive  out-diffusion 
of  Ga  from  the  GaAs  into  the  SiO;.!  J  This  Ga  loss  was  be¬ 
lieved  to  produce  excess  Ga  vacancies  at  the  GaAs  surface. 
Thus,  enhanced  diffusion  of  Si  under  SiO,  encapsulation  is 
similar  to  the  enhanced  diffusion  under  As1'3  and  AsH/ 
overpressures  in  that  both  techniques  produce  excess  Ga  va¬ 
cancies. 

In  another  study,  Kasahara  et  a!.s  showed  that  the 
thickness  of  the  encapsulant  affects  the  diffusion  rate  of  Si  in 
GaAs.  The  diffusion  rate  at  850  °C  under  Si,N4  encapsula¬ 
tion  increased  by  20%  when  the  encapsulant  thickness  was 
increased  from  500  to  1000  A.  Since  the  interfacial  stress 
produced  by  an  encapsulant  scales  directly  with  its  thick¬ 
ness,  this  thickness  dependence  suggested  a  stress-enhanced 
diffusion  rate. 

In  a  similar  study  by  Onuma  et  al.  of  Si  diffusion  under 
Si,N4  encapsulation,  the  measured  diffusion  rate  was  only 
8%  of  the  value  measured  by  Kasahara  et  al:  although  the 
experimental  conditions  were  nearly  identical.  Both  of  these 
measurements  were  made  after  850  °C  anneals  using  a  0.1- 
nm  Si,N4  encapsulant  deposited  by  chemical  vapor  deposi¬ 
tion  techniques.  Since  the  stress  in  Si,N4  films  is  known  to 
depend  on  the  deposition  parameters. ' 1  Si  diffusion  rates  en¬ 
hanced  by  stress  should  show  some  dependence  on  Si,N4 
deposition  conditions.  Different  encapsulant  stresses  could 
explain  the  discrepancy  between  the  results  of  these  two  ex¬ 
periments. 

In  addition  to  the  possibly  different  Si,N4  deposition 
parameters,  the  samples  were  annealed  in  different  atmo¬ 
spheres.  Kasahara  et  al.  measured  the  faster  Si  diffusion  for 


TABLE  I  Summary  of  diffusion  constants  for  Si  in  GaAs  esiracted  from  electrical  measurements  for  a  wide  range  of  post  implantation  annealing  conditions 


Author 

Reference 

Anneal 

ambtem 

Encapsulant 

Anneal 

temperature 

Diffusion 

constant 

Lee  ei  a!. 

4 

H:  over  GaAs/Ga 

None 

8C0'C 

N^ndetectable 

Kasahara  el  at 

5 

5-Torr  AsH,  in  H. 

None 

850  "C 

>  •  10"  cm  /sec 

Kasahara  et  al. 

5 

5-Torr  AsH,  in  H , 

0  1  -ftm  Si,N4 

850 ’C 

.2  .•  10  *  vnv/Nce 

Onuma  et  al 

7 

H; 

0.  f  2-/im  Si  ,Nj 

850  "C 

!  10'  '  cm  - '  sec 

Onuma  et  al 

7 

H 

0  2-/rm  SiO. 

850  C 

2  ■  10  4  cnr/sec 

Kasahara  et  al 

6 

5-Torr  AsH ,  in  H. 

None 

<XXVC 

2  •  10"  4  cm’ 'sec 

Kasahara  et  al. 

6 

0.3 -Torr  AsH,  in  H: 

None 

V00‘C 

!  ■  10"  *  cm  -  /sec 
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>amples  annealed  in  5-mTorr  AsH,  in  H;  while  Onuma's 
!  samples  were  annealed  in  100%  H,.  Previous  results  by  Ka- 
‘  sahara  and  Watanabe"  showed  that  the  St  diffusion  rate  in 
unencapsulated  GaAs  increased  with  AsH ,  partial  pressure. 
Although  the  present  discrepancy  involves  Si,N4  encapsu¬ 
lated  samples,  some  communication  between  the  AsH,  and 
GaAs  may  have  occurred  through  the  encapsulant. 

In  summary,  for  all  the  cases  reviewed  here  where  the 
concentration  of  Si  is  low,  the  diffusion  process  can  be  mo¬ 
deled  using  a  constant  diffusion  coefficient.  It  is  clear  from 
these  experiments,  however,  that  the  value  of  the  Si  diffusion 
constant  in  GaAs  depends  on  the  condition  of  the  GaAs 
surface.  Encapsulants  and  arsenic-bearing  overpressures 
can  vary  the  diffusion  rates  over  an  order  of  magnitude. 
These  annealing  conditions  increase  the  concentration  of  Ga 
vacancies,  suggesting  that  Ga  vacancies  are  an  important 
part  of  the  Si  diffusion  process. 

Results  from  Si  diffusion  experiments  at  high  Si  concen¬ 
trations  exhibited  a  similar  dependence  on  encapsulation 
and  arsenic  overpressures.  However,  the  basic  diffusion  pro¬ 
cess  appears  to  be  different  and  can  be  modeled  with  a  Si  pair 
mechanism.  In  the  next  section,  this  pair  diffusion  mecha¬ 
nism  will  be  developed  fully,  and  solutions  to  the  diffusion 
equation  will  be  presented  that  cover  a  range  of  pairing  con¬ 
ditions. 

III.  A  MODEL  FOR  DIFFUSION  OF  AMPHOTERIC 
DOPANTS  IN  lll-V  SEMICONDUCTORS 

A.  Assumptions 

Recently,  it  was  observed  that  Si  diffusion  into  GaAs 
could  be  modeled  by  a  Si  pair  diffusion  mechanism  when  Si 
was  present  in  large  concentrations^  Here,  we  generalize 
this  pair  mechanism  to  any  amphoteric  dopant  in  a  III-V 
semiconductor  by  solving  the  diffusion  equation  under  a 
wide  range  of  pairing  conditions. 

A  column  IV  dopant  forms  a  donor  when  occupying  a 
column  III  lattice  site  and  forms  an  acceptor  when  occupy¬ 
ing  a  column  V  lattice  site.  Two  column  IV  dopants  residing 
on  nearest-neighbor  donor-acceptor  sites  form  a  neutral 
pair.  It  is  this  neutral  pair  which  is  assumed  to  be  the  mobile 
diffusion  species.  The  concentrations  of  these  dopant  config¬ 
urations  are  written  as  Cd,  C„ ,  and  Cp,  respectively.  The 
total  concentration  of  the  column  IV  dopant  becomes 

CT  =  Cd  +C„  +2C„.  (3) 

Before  proceeding  with  a  diffusion  model  involving  these 
neutral  pairs,  their  concentration  must  by  quantified.  To  cal¬ 
culate  the  percent  of  the  dopant  existing  as  pairs  we  need  to 
make  the  following  assumptions.  First,  neutral  pairs  are  as¬ 
sumed  to  be  in  equilibrium  with  donors  and  acceptors.  This 
is  written  for  Si-doped  GaAs  as 

S*o»  +  SiA,  SiG,  —  SiAs.  W 

The  law  of  mass  action  can  be  applied  to  this  reaction  to 
relate  the  relative  concentrations  of  donors,  acceptors,  and 
pairs.  This  results  in  the  following  pair  equilibrium  equation: 

Cd  Ca=Kp  C„.  (5) 

where  Kp  is  the  pair  equilibrium  constant. 


Second,  if  the  dopant  concentration  is  high  such  that 
the  semiconductor  is  heavily  compensated,  we  can  set  the 
donor  and  acceptor  concentrations  essentialy  equal.  This  is 
written  as 

Cd=C„.  (6) 

This  second  assumption  fails  at  low  doping  levels  where 
most  amphoteric  dopants  preferentially  form  either  donors 
or  acceptors.  The  doping  concentration  where  this  occurs 
depends  upon  the  specific  dopant  and  the  host  lattice  as  well 
as  temperature  and  annealing  conditions. 

Using  Eqs.  (3),  (5),  and  (6)  we  solve  for  the  fraction  of  the 
dopant  that  is  paired,  2CP/CT,  in  terms  of  the  total  dopant 
concentration.  The  result  is 

i£L=i+ifL(l-vTT2  CT/K~).  (7) 

CT  CT 

The  percent  of  the  dopant  paired  is  plotted  in  Fig.  2  versus 
total  dopant  concentration  CT  normalized  to  the  pair  equi¬ 
librium  constant  Kp .  For  CT/Kp  =  1,  27%  of  the  dopant  is 
predicted  to  be  paired  with  the  percent  paired  increasing 
with  total  concentration. 

We  now  wish  to  solve  the  diffusion  equation  assuming 
that  dopant  pairs  dominate  the  diffusion  process.  Under 
these  conditions,  the  dopant  flux  can  be  written  (in  one  di¬ 
mension)  as 

Flux  =  -  2  •  D„  •  ^  .  (8) 

dx 

where  Dp  is  the  pair  diffusion  constant.  Since  the  concentra¬ 
tion  of  pairs  can  be  expressed  in  terms  of  the  total  dopant 
concentration  using  Eq.  (7),  we  can  rewrite  the  flux  equation 
in  a  more  convenient  form  as 


FIG  2.  Percent  of  an  amphoteric  dopant  existing  as  pairs  as  a  function  of 
the  total  dopant  concentration  normalized  to  the  pair  equilibrium  constant. 
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Using  Eq.  |7)  to  evaluate  the  derivative  in  Eq.  1 10)  we  have 

Z>„  =/>,[!-(!  +  2Cr/Kp)-"2].  (1! 

This  effective  diffusion  coefficient  resulting  from  pair  diffu¬ 
sion  is  dependent  upon  the  dopant  concentration.  For  large 
concentrations  [CT  >  Kp ),  the  second  term  in  the  brackets  in 
Eq.  (11)  drops  out.  The  effective  diffusion  coefficient  then 
becomes  constant  and  equal  to  Dp .  As  the  total  dopant  con¬ 
centration  drops,  the  effective  diffusion  coefficient  decreases 
due  to  a  pair  concentration  which  is  dropping  faster  than  the 
total  dopant  concentration. 

B.  Solution  to  the  diffusion  equation 

The  diffusion  equation  can  be  solved  using  the  effective 
diffusion  coefficient  in  Eq.  (11)  and  assuming  a  constant  sur¬ 
face  concentration  Co-  Since  the  value  of  the  pair  equilibrium 
constant  can  vary  over  a  wide  range  depending  on  the  specif¬ 
ic  dopant,  solutions  will  be  presented  for  several  values  of 
Kp .  The  presentation  of  the  solutions  can  be  simplified  by 
noting  that  the  Boltzmann  substitution,  ij  =  x/2 <[7,  can  be 
applied  to  this  diffusion  model.  Thus,  the  depth  scale  is  nor¬ 
malized  to  V/ZT7,  where  Dp  is  the  pair  diffusion  coefficient 
and  t  is  the  diffusion  time.  The  concentration  scale  is  normal¬ 
ized  to  the  surface  concentration  C0.  The  solution  of  the 
diffusion  equation  for  several  pair  equilibrium  constants  is 
shown  in  Fig.  3  with  these  normalized  coordinates. 

For  low  values  of  Kp,  a  deeper  diffusion  is  observed. 
This  is  due  to  the  greater  pairing  which  is  predicted  for  these 
low  values  of  Kp .  To  emphasize  this  point,  we  note  that  for 
Cr  =  O.lCo,  8%  pairing  is  predicted  for  Kp  —  0.5Co  while 
38%  pairing  is  predicted  for  Kp  —  0.05Co.  Thus,  for  the 
same  pair  diffusion  constant  and  diffusion  time,  a  deeper 
diffusion  is  produced  by  the  lower  pair  equilibrium  con¬ 
stants. 

For  all  values  of  Kp ,  Fig.  3  shows  very  steep  diffusion 
fronts.  This  effect  can  be  understood  if  the  diffusion  process 


FIG  3  Solutions  lo  the  diffusion  equation  assuming  a  pair  diffusion  mecha¬ 
nism  and  a  constant  surface  concentration  for  carious  pairing  conditions 


is  visualized  in  the  following  steps.  Some  of  (he  dopant  enters 
(he  lattice  as  neutral  pairs.  These  pairs  are  highly  mobile  and 
quickly  move  to  the  diffusion  front  where  the  total  dopant 
concentration  is  low.  Here  the  pairs  dissociate  into  less  mo¬ 
bile  single  donors  and  acceptors  as  predicted  by  Eq.  (5).  Dif¬ 
fusion  beyond  this  region  then  proceeds  very  slowly  due  to 
the  low  concentration  of  dopant  pairs.  The  result  is  a  very 
steep  diffusion  front.  Deeper  diffusions  produced  by  a  higher 
pairing  percentage  do  not  begin  the  disassociation  process 
until  much  lower  concentrations. 

At  this  time,  two  limitations  of  this  model  should  be 
mentioned.  First,  the  flux  of  single  diffusing  species  may  be¬ 
come  important  at  low  concentrations  where  the  concentra¬ 
tion  of  dopant  pairs  becomes  negligible.  This  model  will  not 
predict  accurately  the  diffusion  at  low  concentrations  since 
mechanisms  not  included  in  the  model  may  become  signifi¬ 
cant.  Second,  the  calculation  of  the  pair  concentration  is  not 
accurate  for  dopant  concentrations  near  the  measured  elec¬ 
tron  concentration.  Here,  the  assumption  of  complete  com¬ 
pensation  breaks  down,  and  the  donor  and  acceptor  concen¬ 
trations  are  no  longer  equal.  The  results  of  the  model  are  not 
affected  significantly  since  correction  of  this  second  approxi¬ 
mation  decreases  the  calculated  pair  concentration  which 
increases  the  slope  of  the  steep  diffusion  front. 

In  summary,  a  diffusion  model  involving  the  formation 
and  rapid  diffusion  of  neutral  dopant  pairs  predicts  concen¬ 
tration-dependent  diffusion  profiles.  The  diffusion  front  is 
very  steep  due  to  the  pair  dissociation  at  low  dopant  concen¬ 
trations.  The  junction  depth  is  dependent  on  the  pair  diffu¬ 
sion  constant,  diffusion  time,  and  pair  equilibrium  constant. 
This  model  is  potentially  applicable  to  all  amphoteric  do¬ 
pants  in  III-V  compounds  where  dopant  pairing  is  likely  to 
dominate  the  diffusion  process. 

IV.  SI  DIFFUSION  IN  GaAs 
A.  Si  pair  equilibrium  constant 

The  pair  diffusion  mechanism  presented  in  the  last  sec¬ 
tion  predicts  diffusion  profiles  that  can  be  shown  to  agree 
with  Si  diffusion  profiles  in  GaAs  at  high  concentrations.' 
Here,  a  pair  equilibrium  constant  of  4x  10'Vcm'  provides 
the  best  fit  for  Si  diffusion  at  1050  *C.  This  gives  a  normal¬ 
ized  pair  equilibrium  constant  of  Kp/C0  =  0.2  for  the  mea¬ 
sured  Si  surface  concentration  of2x  lO’Vcm3.  A  sample  of 
this  fit  is  shown  in  Fig.  4. 

The  application  of  this  pair  mechanism  to  Si  diffusion  in 
GaAs  would  be  supported  strongly  by  an  independent  mea¬ 
surement  of  the  Si  pair  equilibrium  constant.  An  estimate  of 
this  constant  can  be  obtained  from  measurements  of  the  in¬ 
frared  absorption  in  Si-doped  GaAs  due  to  localized  vibra¬ 
tional  modes.  In  a  study  by  Kung  and  Spitzer,13  absorption 
due  to  Si  donors,  acceptors,  and  pairs  was  measured  in  GaAs 
samples  uniformly  doped  with  4x  10"'  Si/cm'.  If  one  as¬ 
sumes  all  Si  atoms  to  have  the  same  absorption  cross  section, 
the  concentrations  of  Si  donors,  acceptors,  and  pairs  can  be 
calculated  from  the  area  of  these  absorption  peaks  and  the 
total  Si  concentration.  These  values  are  found  in  Ref.  12  for 
three  separate  samples  annealed  at  1 100  ’C  and  are  listed  in 
Table  II.  The  pair  association  constant  for  1 100  *C  can  be¬ 
ll 
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FIG.  4.  SIMS  profile  for  a  1030  *C  Si  diffusion  into  GaAs  with  a  calculated 
pair  diffusion  profile  using  Dpt=  1.5  X  10'  lucm;  and  Kr/C0  =  0.2. 


calculated  using  these  values  in  Eq.  (S)  and  is  also  listed  in 
Table  II  for  each  sample.  These  values  are  very  close  to  the 
value  of  Kp  which  gave  the  best  fit  to  the  Si  diffusion  profile 
of  Ref.  3.  The  fact  that  this  independent  measurement  of  the 
pair  equilibrium  constant  agrees  with  the  value  used  to  fit  the 
pair  diffusion  mechanism  to  Si  diffusion  data  strongly  sup¬ 
ports  the  application  of  this  model  to  Si  diffusion  in  GaAs 


B.  Atomic  jump  process  of  SI  pairs 

Mechanisms  for  Si  diffusion  that  involve  motion  of  iso¬ 
lated  silicon  atoms  through  Ga  and  As  vacancies  require 
changes  in  vacancy  type  and  charge  state  that  would  presu¬ 
mably  produce  very  low  diffusion  coefficients.  However,  the 
diffusion  of  a  nearest-neighbor  dopant  pair  is  much  simpler. 
Paired  Si  atoms  can  move  substitutionally  by  exchanging 
sites  with  either  a  Ga  or  As  vacancy.  The  leading  Si  atoms 
moves  into  an  adjacent  vacancy  with  the  trailing  Si  atom 
leaving  behind  the  original  vacancy.  This  process  is  ex¬ 
pressed  by  the  following  equation: 

Sioa  —  SiA,  +  VG,  — *  Va,  +  SiA,  —  SiG,  ■ 

with  a  similar  expression  for  a  pair  exchange  with  an  As 
vacancy. 


Si  pairs  move  through  the  GaAs  lattice  by  alternating 
exchanges  with  Ga  and  As  vacancies.  The  rate  at  which 
these  exchanges  occur  will  depend  upon  the  number  of  Si 
pairs  which  have  the  energy  required  for  the  exchange  and 
on  the  fraction  of  the  time  that  an  adjacent  lattice  site  is 
vacant.  The  first  effect  is  temperature  dependent,  and  the 
second  is  proportional  to  the  Ga  and  As  vacancy  concentra¬ 
tions.  Calculations13  and  experiment14  have  suggested  that 
the  arsenic  vacancy  concentration  exceeds  the  gallium  va¬ 
cancy  concentration  by  at  least  one  order  of  magnitude.  Be¬ 
cause  of  the  lower  Ga  vacancy  concentration,  one  would 
expect  the  waiting  time  for  Ga  vacancies  to  be  much  longer 
than  for  As  vacancies.  This  would  make  the  Si  pair-Ga  va¬ 
cancy  exchange  the  rate-limiting  step  in  the  diffusion  pro¬ 
cess.  Any  moderate  change  in  the  concentration  of  As  va¬ 
cancies  should  not  affect  the  diffusion  rate  while  an  increase 
in  the  Ga  vacancy  concentration  should  result  in  a  corre¬ 
sponding  increase  in  the  overall  diffusion  rate. 

In  an  experiment  by  Vieland1  described  in  Sec.  II,  an 
increase  in  the  arsenic  overpressure  during  diffusion  of  hea¬ 
vily  Si-doped  GaAs  produced  an  increase  in  the  Si  diffusion 
rate.  As  overpressures  are  believed  to  increase  the  concentra¬ 
tion  of  Ga  vacancies.  Thus,  this  result  is  in  agreement  with 
the  suggestion  above  that  Si  pair  exchange  with  a  Ga  va¬ 
cancy  is  the  rate-limiting  step  in  the  pair  diffusion  mecha¬ 
nism.  For  Si  diffusion  into  GaAs  during  rapid  thermal  pro¬ 
cessing,  SiO,  encapsulation  greatly  enhanced  the  diffusion 
rate.  Since  Ga  dissolution  by  SiO.  is  believed  to  produce 
excess  Ga  vacancies,  this  enhancement  also  suggests  that  the 
Si  pair-Ga  vacancy  exchange  is  the  rate-limiting  step  in  the 
Si  diffusion  process.  Thus,  vacancies  are  an  important  consi¬ 
deration  in  the  diffusion  process.  In  the  next  section,  we  con¬ 
sider  the  effect  of  vacancies  on  the  electrical  properties  of  Si- 
doped  GaAs  as  well. 


V.  ELECTRON  COMPENSATION  EFFECTS  IN  Si-DOPED 
GaAs 


A.  Review 


Diffusion  of  Si  into  GaAs  from  a  thin-film  source  intro¬ 
duced  Si  at  levels  well  above  the  measured  electron  concen¬ 
tration.2,3  For  diffusions  at  1050  ’C,  the  electron  concentra¬ 
tion  was  uniform  throughout  the  diffused  region  at 
5  5x  10l8/cm3  while  the  Si  concentration  was  as  high  as 
2x  102O/cmJ  at  the  surface.  This  saturation  effect  is  present 
in  other  doping  processes  as  well.  Davies  ei  a!.'*  reported  a 
maximum  electron  concentration  of  6.5  x  10l8/cm3  for  an- 


TABLE  II.  Relative  infrared  absorption  in  Si-doped  GaAs  due  to  various  Si  configurations  for  three  samples  annealed  at  1 100  *C.  AI»o  listed  are  the  pair 
equilibrium  constants  calculated  with  Eq.  Ul. 


Sample 

Relative  absorption 

Pair  equilibrium 
constani  K. 

Si,.. 

Si,. 

Sin. -Si  „ 

4 

212 

96 

45 

4.5  X  lO'Vcm’ 

5 

211 

82 

45 

4  Ox  10‘Vcm 

6 

225 

91 

42 

4.9  x  10'Vcm ' 

12 


\ 

nealed  Si  implants  where  the  peak  Si  concentration  was 
2  X  lO'Vcm  Chat  tt  at.  showed  that  electron  concentra¬ 
tions  in  St-doped  molecular  beam  epitaxy  GaAs  peak  at 
6X  lO'Vcm’  while  a  group  from  Rockwell17  obtained  a 
higher  activation  by  changing  the  growth  conditions. 

These  electron  saturation  effects  have  been  explained 
by  a  number  of  methods.  There  are  electronic  degeneracy 
effects  in  semiconductors  that  predict  some  form  of  free- 
carrier  saturation  at  high  doping  levels.  Effects  due  to  the 
chemical  solubility  of  impurities  also  have  been  used  to  ex¬ 
plain  free-carrier  saturation  effects  in  GaAs.  In  Te-doped 
GaAs,  the  formation  of  a  stable  precipitate.  Ga;Te„  is  be¬ 
lieved  to  prevent  the  incorporation  of  single  substitutional 
Te  above  levels  of  lO'Vcm3.18  For  doping  by  Si,  Ge.  and  Sn. 
free-electron  saturation  can  be  attributed  to  dopant  precipi¬ 
tation.  However,  since  these  dopants  are  amphoteric,  the 
formation  of  compensating  acceptors  and  neutral  donor-ac¬ 
ceptor  pairs  can  produce  the  free-electron  saturation  ob¬ 
served  under  heavy  doping  conditions.  Here,  we  propose  a 
mechanism  that  predicts  the  formation  of  compensating  ac¬ 
ceptors  and  pairs  at  high  electron  concentrations  w  hich  may 
be  a  significant  part  of  the  electron  saturation  effect  found  in 
heavily  Si-doped  GaAs.  The  mechanism  is  applicable  to  oth¬ 
er  amphoteric  dopants  in  III-V  compounds  as  well. 

B.  Compensation  mechanism  and  predicted  saturation 
profiles 

The  equilibrium  between  single  Si  and  vacancies  can  be 
expressed  by  the  following  reaction: 

Sio.  +  y:>  +  2e“  s=*  SiAs  +  V <r,  -t  2h  ~  1 12) 

Based  on  this  compensation  reaction,  the  site  location  of  Si  in 
GaAs  will  be  affected  by  the  electron  and  hole  concentra¬ 
tions  as  well  as  As  and  Ga  vacancy  concentrations.  To  model 
the  compensation  in  Eq.  (12),  all  donors  |Si0a  and  V ’Ai)  and 
acceptors  (SiA,  and  VQt )  are  treated  as  completely  ionized. 
In  view  of  the  Fermi-level  position  and  the  temperatures 
used  in  this  study,  this  is  a  good  assumption  for  all  species 
except  VAi  which  does  not  have  a  well-characterized  donor 
level.  Whether  the  As  vacancies  are  neutral  or  ionized  will 
not  change  the  basic  results  of  the  compensation  mecha¬ 
nism. 

If  equilibrium  exists  among  the  species  in  Eq.  1 12),  the 
law  of  mass  action  can  be  used  to  calculate  their  relative 
concentrations.  The  result  is 

[SiG+.  ]  •  [  V:>  ]n-  =  Kc-  [SiA5  ]  ■  [  FG-a  ]  -p-.  (13) 

where  the  brackets  [  ]  denote  concentration.  The  equilibri¬ 
um  constant  Kc  is  related  to  the  activation  energy  for  the 
reaction  in  Eq.  1 121.  Free-carrier  saturation  is  inherent  in  this 
mechanism  and  can  be  seen  by  observing  the  effect  of  chang¬ 
ing  n  or p  on  the  distribution  of  Si  between  Ga  and  As  sites. 
We  assume  for  now  that  the  vacancy  concentrations  are 
fixed.  By  increasing  n  in  Eq.  (13),  Si  shifts  from  Ga  sites  to  As 
sites  to  maintain  equilibrium.  This  reduces  the  donor  con¬ 
centration  and  increases  the  acceptor  concentration,  de¬ 
creasing  n.  The  reverse  effect  occurs  when  p  is  increased  in 
Eq.  (13). 

The  sensitivity  of  the  compensation  to  the  electron  con¬ 


i'avcr  y 

cent  ration  can  be  extracted  from  Eq. :  1 3i  by  solving  tor  the  si 
acceptor  to  donor  ratio  y.  This  is  also  called  the  compensa¬ 
tion  ratio  and  is  shown  below: 


[SiAl!  =  _L_  [  V  C.  ]  fn_V 

[Sio,]  KAT)  [KGa]  \nj 


1 14) 


Thus,  the  compensation  ratio  is  predicted  to  be  proportional 
to  the  fourth  power  of  the  electron  concentration.  These  re¬ 
sults  are  based  on  positively  charged  As  vacancies.  The  com¬ 
pensation  ratio  will  be  proportional  to  the  third  power  of  the 
electron  concentration  if  the  As  vacancies  are  neutral. 

In  some  liquid  phase  epitaxy  layers.  Si  doping  produces 
p-type  conductivity.  Under  these  conditions 
[SiT,  ]  >  [Sic,  ]  and  y  >  1.  Growth  occurs  under  Ga  melts 
which  produce  larger  As  to  Ga  vacancy  concentration  ratios 
than  other  epitaxial  processes.  According  to  Eq.  i  1 3),  these 
conditions  increase  the  compensation  ratio.  This p-type  con¬ 
ductivity  for  Si  doping  is  rare,  however.  For  most  epitaxial 
growth  and  annealing  conditions  using  low  Si  concentra¬ 
tions,  only  Si  donors  are  produced  resulting  in  nearly  un¬ 
compensated  .n-type  conduction.  For  these  conditions, 
y  <  1 .  For  larger  Si  concentrations,  conduction  is  still  n  type 
but  with  higher  compensation. 

From  Eq.  (13),  the  compensation  ratio  is  predicted  to 
increase  with  the  fourth  power  of  the  electron  concentration 
to  the  intrinsic  concentration  ratio,  ( n/n ,  |4.  For  Si  concen¬ 
trations  and  anneal  temperatures  where  the  electronic  con¬ 
duction  is  intrinsic  (n  =  n, ),  the  compensation  ratio  y  is  pre¬ 
dicted  to  be  independent  of  donor  and  acceptor 
concentrations.  Under  these  intrinsic  conditions,  the  Si  and 
donor  concentrations  are  usually  found  to  be  equal,  imply¬ 
ing  y  <  1.  However,  as  the  doping  increases  and  semicon¬ 
ductor  becomes  extrinsic  (n  >  n, ),  the  compensation  ratio  is 
predicted  to  be  increased  by  the  [n/n,  I4  term.  At  this  point, 
additional  Si  doping  produces  acceptors  as  well  as  donors. 
As  the  compensation  ratio  approaches  unity  (y— -  l),  addi¬ 
tional  Si  doping  is  predicted  to  incorporate  equally  as  donors 
and  acceptors  effectively  limiting  the  electron  concentra¬ 
tion. 


Defining  n  =  n when  y  —  1  in  Eq.  (14)  and  solving  for 
,  gives 


»4.,=«<-(/fc(7').[F0-a]/[F*]|1'4.  (15) 


Thus,  the  saturated  electron  concentration  can  be  shown  to 
depend  upon  terms  which  are  temperature  dependent.  In 
addition,  this  saturation  value  increases  with  the  fourth  root 
of  the  ratio  of  vacancy  concentrations.  At  a  fixed  tempera¬ 
ture  it  is  possible  to  vary  this  ratio  by  changing  annealing  or 
epitaxial  growth  conditions  and  thus  influence  the  saturated 
electron  concentration.  However,  the  level  of  electron  satu¬ 
ration  is  predicted  to  be  weakly  dependent  on  these  condi¬ 
tions.  For  example,  the  ratio  of  Ga  to  As  vacancy  concentra¬ 
tions  must  be  changed  by  a  factor  of  16  to  produce  a 
tactor-of-2  change  in  the  saturated  electron  concentration  of 
Eq.  115). 

The  shape  of  this  saturation  versus  Si  concentration  can 
be  determined  by  solving  for  the  room-temperature  electron 
concentration.  n„  =  [ Si,;,  ]  —  [SiA>  ].  This  differs  from  n , 
the  electron  concentration  at  the  annealing  temperature,  by 
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FIG.  5.  Room-lemperaiureeJectron  concentration  vs  roiaJ  dopant  concen¬ 
tration  calculated  from  Eq.  1 16)  for  various  values  of  K,  ■ 


the  hole  concentration,  p  =  nf/n.  The  hole  concentration  at 
room  temperature  can  be  ignored  for  most  doping  levels. 
The  following  relations  between  the  electron  concentration 
and  the  various  Si  species  will  be  used: 


=  [Sio.]  -  [SiA.  ],  (16a) 

«+  [Si*,]  =/>+  [Sio.  ],  (16b) 

l Si]  =  [Sio,]  +  [SiAs]  +2  •  [SiG,-SiAs],  (16c) 
[Sio, )  •  [SiA4  ]  =  Kp  -  [SiG>  —  SiA,  ] ,  ( 16dl 

««/«/  =  (*,  •  [  I'o.  ]/[  ]),/4  =  Kv.  (16e| 


Equations  (16c)  and  ( 1 6d)  include  the  contribution  of  Si  near¬ 
est-neighbor  pairs  in  the  compensation  process.  Similar 
equations  were  used  in  Sec.  Ill  for  the  pair  diffusion  model. 
We  will  solve  this  set  of  equations  for  rtn  vs  [Si]  using  differ¬ 
ent  values  of  AT„  in  Eq.  (16e).  Note  that  we  have  ignored  the 
contribution  of  the  vacancies  in  the  charge  neutrality  equa¬ 
tion  (16b).  Vacancy  concentrations  are  taken  to  be  in  the 
lO'Vcrn3  range  or  less  and  should  be  considered  for  Si  con¬ 
centrations  in  this  range. 

The  solutions  to  Eqs.  ( 16aH  1  he)  are  plotted  in  Fig.  5  for 
K,  =5,10,20  with  an  intrinsic  carrier  concentration  of 
n,  =  I  x  10'VcmJ  corresponding  to  1000  ”C.  In  each  curve, 
there  is  a  sharp  saturation  at  a  level  determined  by  A',  with 
negligible  compensation  for  lower  Si  concentrations.  From 
Eq.  (14),  we  see  that  the  compensation  ratio  for  these  low  Si 
concentrations  is  equal  to  A'," 4.  which  is  nearly  zero  for  the 
values  of  K ,  in  Fig.  5. 

We  have  selected  values  of  A',,  that  produce  saturation 
curves  typical  of  Si  doping  in  GaAs.  Thedirect  calculation  of 


A'„  as  a  function  of  temperature  and  annealing  conditions  is 
beyond  the  scope  of  this  paper.  However,  the  solutions  in 
Fig.  5  do  suggest  that  the  autocompensation  of  amphoteric 
dopants,  driven  by  an  increasing  electron  concentration,  can 
account  for  the  saturated  free-eiectron  concentrations  wide¬ 
ly  observed  in  GaAs.:  V15'17 

VI.  SUMMARY 

The  amphoteric  nature  of  Si  in  GaAs  has  been  used  to 
develop  diffusion  and  electrical  compensation  mechanisms. 
The  diffusion  mechanism  is  based  on  nearest-neighbor  do¬ 
nor-acceptor  pairs.  General  solutions  were  presented  that 
predicted  abrupt  diffusion  fronts  for  a  wide  range  of  pairing 
conditions.  Experiments  support  the  application  of  this 
mechanism  to  Si  diffusion  in  GaAs  at  high  concentrations. 

A  compensation  mechanism  for  amphoteric  dopants 
was  developed  as  well.  The  compensation  process  is  driven 
primarily  by  the  free-electron  concentration.  Nearly  com¬ 
plete  compensation  is  predicted  for  large  dopant  concentra¬ 
tions. 
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ABSTRACT 

Epitaxial  regrowth  is  investigated  for  layers  of  InSb  and  GaAs 
amorphized  by  liquid  nitrogen  temperature  ion  implants.  Experimental 
criterion  for  amorpnization  and  limitation  of  regrowth  are  correlated 
with  damage  calculated  by  a  Boltzmann  transport  equation  approach  to  ion 
implant  modeling.  Conditions  for  complete  epitaxial  regrowth  as  deter¬ 
mined  by  channeled  Rutherford  backscattering  spectrometry  are  presented 
with  evaluation  of  residual  defects  by  transmission  electron  microscopy. 
Electrical  activation  following  regrowth  of  GaAs  is  reported  and  corre¬ 
lated  with  calculated  damage  profiles. 


INTRODUCTION 

An  understanaing  of  defects  resulting  from  ion  implant  damage  in 
compound  semiconductors  is  vital  to  employing  ion  implantation  effectively 
for  ir  detector  arrays  and  high-speed,  high-density  III-v  integrated 
circuits.  A  technique  has  been  developed  by  Christel  [1]  for  the  modeling 
of  lattice  displacements  and  stoichiometric  imbalances  resulting  from 
implant  damage  in  compound  semiconductors.  This  technique,  the  Boltzmann 
transport  equation  .BTE)  approach  to  ion  implantation  in  multilayered 
structures  and  compound  semiconductors,  is  used  in  this  work  to  correlate 
experimental  observations  with  physical  damage  mechanisms.  Damage  pro¬ 
files  from  Rutherford  Dackscatteri ng  spectrometry  (RBS)  and  transmission 
electron  microscopy  iTE«)  are  correlated  with  features  of  BTE  calculations 
to  develop  a  physical  understanding  of  the  defect  formation  processes  and 
allow  for  prediction  of  results  over  a  spectrum  of  implant  parameters . 

The  particular  application  of  correlation  of  experiment  with  BTE 
calculations  addressed  in  this  work  is  the  amorphi zation  and  solid  phase 
epitaxial  regrowtn  of  compound  semiconductors.  The  motivation  for  an 
amorphizing  ion  ’-plant  followed  by  solid  phase  epitaxy  (SPE)  is  that  it 
nas  been  known  fo'  some  time  [2]  that  in  silicon  this  technique  results 
in  activation  of  'mplanted  dopants  at  significantly  lower  temperatures 
(<6Q0°C)  than  those  required  to  remove  implant  defects  when  the  target  is 
not  amorphized.  *"e  lower  temperature  required  for  the  anneal  and 
activation  of  an  amorphizing  implant  would  be  attractive  for  compound 
semiconductor  devices  from  the  standpoint  of  minimizing  material  decom¬ 
position  as  well  as  reducing  impurity  redistribution,  which  degrades 
the  definition  of  cev ice  geometries.  SPE  of  GaAs  following  amorphization 
was  demonstrated  :v  Williams  and  Austin  [3]  and  has  recently  been 
investigated  by  ethers  [4-7]  but  electrical  activation  has  not  been 
reported.  In  the  -'ol  lowing  sections  amorphization  criteria  are  proposed 
for  InSb  and  GaAs  vol  lowed  by  observations  or.  the  correlation  of  SPE 
results  with  BTE  calculations. 


AMORPHIZATION 

Two  models  for  the  threshold  of  amorphization  in  column  IV  semicon¬ 
ductors  have  been  proposed.  These  are  the  damage  densi‘.y  model  of  Vook 
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[8]  and  the  lattice  displacement  model  of  Swanson  et  al.  [9].  The  damage 
density  criterion  proposes  that  the  energy  deposited  in  nuclear  stopping 
processes  must  exceed  a  threshold  value  for  amorphization  to  occur.  For 
silicon,  at  sufficiently  low  temperature  to  avoid  in-situ  annealing,  the 
threshold  has  been  determined  to  be  6  to  10xlo20((eV/cm^  [8].  The  lattice 
displacement  criterion  proposes  that  the  lattice  will  relax  to  an  amor¬ 
phous  phase  when  a  specific  fraction  of  the  lattice  has  been  displaced. 
Swanson  [9]  has  calculated  a  lattice  displacement  threshold  of  2  to  4t 
for  Ge  and  Christel  [10]  has  determined  a  threshold  of  101  for  Si  by 
correlation  of  8TE  calculations  with  experiment. 

Profiles  of  damage  density  and  lattice  displacement  as  calculated  by 
the  BTE  approach  have  been  compared  with  the  results  of  amorphizing  im¬ 
plants  into  InSb  and  GaAs.  An  example  of  this  correlation  is  shown  in 
Fig.  1  for  a  280KeV  Cd  implant  into  InSb.  Channeled  RBS  spectra  were  used 
to  determine  amorphous  layer  thickness.  For  these  measurements  2.2MeV 
He*  particles  were  channeled  in  the  (111)  JnSb  and  (100)  GaAs.  The 
damaged  layer  is  indicated  by  a  high  yield  of  backscattered  He*.  The  He* 
were  detected  at  an  angle  of  5  degrees  from  the  surface  for  shallow 
implants  to  enhance  depth  resolution.  Comparing  similar  calculations  for 
Cd  implants  into  InSb  at  energies  ranging  from  100  to  280KeV  and  doses  of 
lO*^  to  Ifl*^  ions/cm?,  the  damaging  energy  density  necessary  for  complete 
ancrphizat  ion  of  InSb  was  found  to  lie  in  the  range  of  2  to  3x10^0  keV/ 
cm2,  the  fact  that  this  value  is  less  than  that  observed  for  Si  is  consis¬ 
tent  with  the  cohesive  energy  of  InSb  bonds  ( 1  . 4e V )  being  less  than 
cohesive  bond  energy  of  Si  (2.3eV)  [II],  The  damage  density  threshold 
for  GaAs  has  been  determined  to  oe  on  the  order  or  1.5  to  1 lxlo20keV/cm2 
for  Si  implants  at  50  to  ISOKeV.  The  cohesive  bond  energy  of  GaAs  is  re¬ 
ported  to  be  1.6eV  [11].  For  the  lattice  displacement  calculation  a  va  ue 
must  be  assumed  for  the  energy  necessary  to  displace  an  atom  from  its 
lattice  site  which  would  correspond  to  four  cohesive  bond  energies  [10]. 


CEi’TH  til 

Figure  1  Correlation  of  anorpnous 
layer  thickness  measured  Cy  98S 
with  damage  density  computed  by 
the  BTE  approach  for  a  230KeV/ 
3x10*2  Cd  implant  into  lnSB. 


Figure  2  Correlation  of  BBS  of 
InSb  damaged  layer  as  implanted 
and  after  anneal  with  BTE  calcu¬ 
lations  of  damage  profiles. 
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Correlation  of  BTE  calculations  with  RBS  measurements  of  amorphous  layer 
thickness  results  in  a  lattice  displacement  criterion  of  2  to  31  for  Insb 
and  3  to  164  for  GaAs  using  S.6eV  and  6.4eV,  respectively,  for  displace¬ 
ment  energies. 

Both  the  damage  density  and  lattice  displacement  thresholds  deter¬ 
mined  for  GaAs  vary  By  nearly  an  order  of  magnitude.  This  is  believed  to 
be  the  result  of  annealing  occurring  during  implantation  though  all 
implants  were  done  at  liquid  nitrogen  temperature  in  order  to  minimize 
this  effect.  As  discussed  by  Vook  and  stein  [12],  the  production  rate  of 
defects  in  the  approach  to  the  amorphous  transition  is  the  result  of  a 
balance  between  the  introduction  By  the  ion  flu*  and  the  removal  by 
thermal  annealing  processes.  A  complicating  consideration  is  that  while 
an  increase  in  ion  dose  rate  will  increase  the  defect  introduction  rate 
linearly,  it  will  also  cause  target  heating  which  can  have  an  exponential 
effect  on  the  defect  removal  rate. 


SOLID  PHASE  EPITAXY 

After  rendering  an  InSb  or  GaAs  surface  layer  amorphous,  SPE  was 
performed  to  recover  the  crystal  structure.  RBS  spectra  for  InSb  before 
and  after  regrowtn  are  shown  in  fig.  2  below  calculated  damage  profiles. 
In  this  case  regrowth  is  not  complete,  stopping  near  the  depth  at  which 
a  transition  occurs  from  net  excess  In  and  SB  to  net  vacancies.  For  Cd 
implants  at  100,  220  and  280  Kev  regrowth  appeared  to  start  rapidly  (>1 
Angstrom/sec)  and  then  stop  (<.01  Angstrom/sec)  after  30s  at  280C.  Corre¬ 
lating  the  depth  at  which  regrowth  stops  with  BTE  calculations,  indicates 
that  this  occurs  at  a  net  vacancy  concentration  on  the  order  of  lO^cm’3. 
A  stopping  of  the  regrowth  due  to  vacancies  indicates  a  limi'ation  of  dif¬ 
fusion  of  In  and  So  from  the  depths  to  which  they  are  recoiled,  out  to  the 
near  surface  region  of  net  vacancies.  This  limitation  in  the  rate  of  SPE 
could  permit  time  for  nucleation  of  crystallites  at  the  surface  or  within 
the  damaged  layer.  This  could  stop  SPE  short  of  reaching  the  surface  as 
shown  by  Williams  et  al.  [13],  Random  orientation  of  these  crystallites 
would  cause  the  damaged  layer  to  regrow  in  a  polycrystalline  structure. 
The  slope  of  the  deep  damage  edge  in  the  RBS  data  of  Fig.  2  would  support 
the  argument  of  crystallite  formation.  The  slope  of  the  edge  before 
anneal  is  consistent  with  an  abrupt  amorphous  edge  convolved  with  the  sys¬ 
tem  resolution  ana  energy  straggle.  In  the  post -anneal  case,  however,  the 
slope  indicates  fiat  the  transition  is  no  longer  abrupt  but  graded  due  to 
the  nucleation  of  misoriented  crystallites. 

Regrowth  experiments  on  GaAs  also  demonstrated  a  nonlinear  growth 
rate  characterized  By  initial  rapid  growth  that  slows  and  then  stops  at  a 
given  anneal  temperature.  Figure  3  shows  RBS  data  for  GaAs  as  amorphized 
with  Si  at  SOKeV  and  regrown  at  550C.  In  addition  to  regrowth  stopping 
after  10  minutes  at  550C,  arsenic  loss  near  the  surface  can  be  observed 
as  indicated  by  the  lower  RBS  yield  at  the  surface  in  the  post-anneal 
channeled  spectrum  and  random  spectrum  (not  shown).  To  prevent  arsenic 
loss  during  the  anneal  a  thin  film  of  aluminum  was  evaporated  onto  the 
surface  of  a  sample  that  was  then  processed  with  the  same  implant  and 
anneal  conditions.  The  aluminum  encapsulated  sampled  regrew  completely 
as  indicated  by  cne  MBS  data  of  Fig.  3  showing  the  post-anneal  spectrum 
coincident  with  tne  spectrum  for  unimplanted  material.  Plan  view  bright 
field  TEM  and  electron  diffraction  patterns  are  shown  for  the  aluminum 
capped  sample  as-mplanted  and  after  anneal  in  Fig.  4(a)  and  4(b),  respec¬ 
tively.  In  the  as-implanted  case  amorphous  rings  are  observed  in  the 
diffraction  pattern  superimposed  on  the  pattern  for  the  (100)  substrate. 


17 


Paper  #2 


t 


Figure  3  Channeled 
BBS  spectra  for  GaAs 
amorphized  by  50KeV 
Si  and  regrown  at 
SSOC  for  10  min. 
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After  the  anneal,  residual  dislocation  loops  are  apparent.  BTE  calcula¬ 
tions  of  recoil  and  damage  profiles  for  the  aluminum  capped  sample  are 
shown  in  Fig.  5.  It  may  be  noted  that  this  condition  in  which  complete 
regrowtn  was  achieved  corresponds  to  no  net  vacancies  near  the  surface 
due  to  aluminum  recoils. 


A1  j  GaAs 
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Figure  A  Plan-view  bright-field  TEM 
and  electron  diffraction  patterns  for 
GaAs  (a)  as  implanted  and  (Pi  after 
550C  for  10  min. 


Figure  i  Correlation  of  elec 
trical  activation  pro'iles  of 
regrown  GaAs  with  BTE  calcula 
tions  of  recoil  profiles 
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Electrical  activation  analysis  By  van  der  Pauw-HaU  measurements 
were  performed  on  an  aluminum  capped  sample  of  (100)  undoped,  LEC  GaAs 
and  the  results  are  shown  in  Fig.  5.  A  low  level  of  n-type  electrical 
activation  was  observed  to  the  depth  of  the  initial  amorphous  layer  in¬ 
dicating  that  some  of  the  implanted  Si  replaced  Ga  atoms  in  the  low  tem¬ 
perature  regrowth  process.  The  level  of  activation  near  the  surface  is 
suppressed,  possibly  due  to  A1  preferentially  growing  onto  Ga  sites  as 
opposed  to  Si  in  this  region  where  A!  concentration  is  high.  The  high 
mobility  observed  is  consistent  with  a  low  concentration  of  residual 
electrically  active  defects. 


CONCLUSION 

Amorphi  rat  ion  hy  ion  implantation  and  subsequent  SPE  of  InSb  is 
reported  for  the  first  time.  The  threshold  for  amorphi zat i on  of  InSb  ana 
GaAs  in  terms  of  damage  deposited  in  nuclear  stopping  processes  has  been 
determined  by  correlating  RBS  measurements  to  BTE  calculations  of  damage. 
A  limitation  of  SPE  in  GaAs  and  InSb  appears  to  be  correlated  with  the 
transition  from  a  net  excess  deep  in  the  implanted  material  to  net  vacan¬ 
cies  near  the  surface.  This  atomic  imbalance  is  due  to  recoils.  GaAs 
layers  under  an  A!  cap  have  been  amorphized  and  regrown  completely.  In 
tnis  case,  A1  recoiled  into  the  GaAs  compensates  for  Ga  and  As  displaced 
near  the  Surface  to  result  in  no  net  atomic  vacancies.  A  low  level  of  n- 
type  electrical  activation  was  observed  in  the  regrown  GaAs  with  mobility 
indicative  of  few  residual  electrically  active  defects.  With  the  implant 
damaye  removed  by  the  low  temperature  regrowth  a  rapid  anneal  couid  be 
employed  for  further  Si  activation. 
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Extended  defects  in  GaAs  are  investigated  following  epitaxial  regrowth  of  amorphous  layers. 

GaAs  surface  layers  were  amorphized  by  Si*  implants  at  liquid-nitrogen  temperature.  Anneals 
were  performed  for  4  s  to  20  min  from  1 50  to  885  'C.  Rutherford  backscattering  spectrometry 
and  transmission  electron  microscopy  were  used  to  evaluate  the  results  of  annealing.  Complete 
solid-state  epitaxy  occurs  rapidly  at  low  temperature.  Stacking  faults  and  microtwins 
surrounded  by  dislocation  networks  extend  to  the  surface  following  regrowth.  The  dislocations 
anneal  at  varying  rates  over  intermediate  temperature  ranges  (200-700  *C),  and  the 
microtwins  climb  out  at  higher  temperatures  (  >  700  ‘C).  Defect  depth  profiles  are  correlated 
with  damage  and  stoichiometric  imbalances  computed  by  the  Boltzmann  transport  equation 
approach  to  ion-implant  modeling. 


It  is  well  known  that  solid  phase  epitaxy  (SPE)  at  a 
relatively  low  temperature  (  -600*C  for  a  Si  crystal)  per¬ 
mits  ion-implanted  dopants  to  be  incorporated  into  substitu¬ 
tional  sites  in  a  Si  crystal.  This  process  was  first  demonstrat¬ 
ed  for  silicon  by  Crowder  in  1970.'  For  the  compound 
semiconductor  GaAs  a  low-temperature  process  for  impuri¬ 
ty  incorporation  is  attractive  because  it  provides  both  re¬ 
duced  material  decomposmon  and  reduced  impurity  redis¬ 
tribution  during  annealing.  SPE  of  GaAs  was  demonstrated 
by  Williams  and  Austin  in  1980:  and  has  recently  been  inves¬ 
tigated  by  others  ''7  Although  SPE  in  GaAs  has  been  dem¬ 
onstrated  at  very  low  temperatures7  (180*0.  extended  de¬ 
fects  are  often  observed  in  the  implanted  layer  following 
SPE  In  this  letter,  Rutherford  backscattering  spectrometry 
(RBS)  and  high  resolution  transmission  electron  micros¬ 
copy  ( HRTEM )  results  for  Si-implanted  GaAs  are  present¬ 
ed  with  a  correlation  to  lattice  damage  calculated  by  the 
Boltzmann  transport  equation  (BTE)  approach  to  ion  im¬ 
plantation  modeling."'1"  Based  on  this  correlation  a  source 
of  microtwin  nucleation  is  presented,  extrinsic  dislocation 
loops  are  tentatively  identified,  and  an  amorphization 
threshold  is  derived. 

Silicon  ions  were  implanted  into  undoped  semi-insulat¬ 
ing  liquid  encapsulated  Czochralski  ( 100)  GaAs  to  a  dose  of 
10''cm  7  at  100  keV.  The  GaAs  target  was  cooled  to  liquid- 
nitrogen  temperature  during  implantation  to  minimize  in 
situ  annealing.  After  rendering  a  GaAs  surface  layer  amor-’ 
phous,  SPE  was  performed  to  recover  the  crystal  structure. 
Anneals  were  conducted  over  a  temperature  range  of  ISO- 
885  *C  in  forming  gas.  Anneal  times  ranged  from  4  s  to  30 
min.  For  HRTEM  analysis  cross-sectional  specimens  were 
prepared  by  a  technique  similar  to  that  described  by  Brav¬ 
man"  with  modificat  ons  to  minimize  specimen  heating 
which  could  lead  to  unwanted  annealing. 

Depth  profiles  of  the  BTE  calculations  for  displaced  Ga 
or  As  atoms,  shown  at  the  lop  of  Fig.  1.  are  computed  to  be 


more  than  an  order  of  magnitude  greater  than  the  Si  concen¬ 
tration  over  much  of  the  profile.  Although  computed  pro¬ 
files  for  Ga  and  As  differ  slightly  due  to  their  difference  in 
mass,  they  are  shown  in  Fig.  1  as  a  single  curve  since  differ¬ 
ences  are  generally  less  than  a  few  percent.  The  profile  of  Ga 
or  As  displaced  was  calculated  assuming  a  10-eV  displace¬ 
ment  energy.  Though  the  appropriate  value  for  displace¬ 
ment  energy  is  subject  to  speculation,  use  of  5  and  20  eV  was 
found  to  result  in  only  a  20%  change  in  the  displacement 
concentration  profile.  The  recoiling  of  Ga  and  As  from  the 
near-surface  region  results  in  a  net  vacancy  concentration 
there  while  Ga  and  As  atoms  pile  up  deeper  into  the  sub¬ 
strate  producing  a  local  net  excess.  Due  to  the  difference  in 
Ga  and  As  masses,  the  Ga  atoms  are  recoiled  further,  result¬ 
ing  in  the  surface  region  being  slightly  As  rich  while  the 
deeper  region  is  slightly  Ga  rich. 

Transmission  electron  microscopy  (TEM)  on  the  As- 
implanted  material  shown  in  Fig.  2  reveals  an  amorphous 
region  extending  from  the  surface  to  a  depth  of  950  A.  The 
amorphous  region  was  indicated  by  an  absence  of  thickness 
contours  on  tilting  the  specimen  and  by  the  rings  in  the  se¬ 
lected  area  diffraction  pattern.  From  950  to  1250  A  a  heavily 
defected  transition  region  is  observed  to  occur  between  the 
amorphous  layer  and  the  single  crystal  bulk.  The  threshold 
for  amorphization  was  determined  in  terms  of  two  BTE  cri¬ 
teria.  Comparing  the  depth  of  the  transition  region  to  the 
computed  displacement  profile  indicates  that  the  fractional 
lattice  displacement  necessary  for  a  crystalline  to  amor¬ 
phous  transition  in  GaAs  is  in  the  range  of  1 8-25%.  A  simi¬ 
lar  correlation  was  made  with  a  computed  profile  of  the  ener¬ 
gy  deposited  in  nuclear  stopping  processes,  indicating  an 
amorphization  threshold  in  terms  of  this  criterion  of  1.1- 
1.7x  1071  keV/cm'.  This  corresponds  to  50-77 eV  per  GaAs 
molecule  for  Si  implantation,  significantly  higher  than  the 
11  +  2  eV  reported  by  Sands  et  at.  for  Se  implantation. 

Following  SPE,  the  near-surface  region  is  seen  to  be 
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FIG  l  Depth  profiles  of  Bollimann  transport  equation  calculations  for  a 
10'*  Si/cm:  implant  at  100  KeV  into  GaAs 

highly  defected  in  the  cross-sectional  TEM  of  Fig.  2.  The 
selected  area  diffraction  pattern  indicates  the  regrowth  was 
complete  but  a  high  density  of  stacking  faults  and  mi¬ 
crotwins  remain.  Stacking  faults  were  indicated  by  the 
streaks  in  the  diffraction  pattern  and  microtwins  were  indi¬ 
cated  by  the  faint  diffracted  beams  displaced  one-third  of  the 
distance  between  primary  diffracted  beams. 

A  lattice  image  for  this  specimen  is  shown  in  Fig.  3 
viewed  in  a  <  1 10)  direction.  For  this  image  four  1 1 1  and  two 
200  reflections  were  used  in  conjunction  with  the  undiffract- 
ed  beam.  Here  microtwins  are  clearly  evident  with  their 
characteristic  boundaries  on  (111)  facets  intersecting  the 
( 100)  material  surface  at  an  angle  of  54.7°  The  HRTEM 
demonstrates  that  crystal  recovery  with  low  defect  densities 
occurs  in  the  transition  region.  In  regrowth  of  the  amor¬ 
phous  region,  however,  stacking  faults  and  microtwins  nu¬ 
cleate  and  propagate  to  the  surface. 

Although  the  average  net  Ga  and  As  concentrations  are 
too  low  to  account  for  sufficient  antisite  defects  to  nucleate 
all  the  stacking  faults  observed,  local  variations  in  the  stoi¬ 
chiometry  around  damage  clusters  could  account  for  this. 
BTE  calculations  indicate  that  greater  than  99%  of  the  dis¬ 
placed  atoms  at  the  amorphous  to  crystalline  interface  recoil 
less  than  5  A.  The  density  of  stacking  faults  and  microtwins 
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FIG  2.  Cross-sectional  TEM  of  GaAs  (a)  as  implanted  with  lOOkeVSi  to  a 
dose  of  I01'  cm" ’and  after  anneals  of  (b)  400  °C  for  10  s.  (c)  600  "C  for  30 
s,  (d  I  700  "C  for  30  s,  and  (e)  800  *C  for  30  s. 


in  the  regrown  region  does  correlate  with  the  high  density  of 
total  displacements.  A  high  density  of  residual  dislocation 
loops  is  observed  in  TEM  of  the  completely  regrown  materi¬ 
al.  The  depth  at  which  these  are  observed  correlates  with  the 
Si  peak  and  the  net  Ga  and  As  excess  in  the  BTE  calcula¬ 
tions,  indicating  that  these  may  be  extrinsic  loops  formed  by 
Si,  Ga,  and  As  precipitation. 

Channeled  RBS  spectra  are  shown  in  Fig.  4  for  anneals 
at  600,  700,  and  800  °C.  These  anneals  were  intended  to  re¬ 
move  residual  defects  following  SPE  at  400  °C  for  10  s.  The 
depth  scale  shown,  though  based  on  a  surface  approximation 
of  electronic  stopping  power,  correlates  well  with  cross-sec¬ 
tional  TEM.  An  anneal  at  600  °C  for  4  s  resulted  in  both  a 
reduction  in  the  depth  to  which  the  damage  extends  and  a 
reduction  in  the  height  of  the  near-surface  yield.  Anneals  for 
longer  times  at  600  °C  can  be  seen  to  further  reduce  the  depth 
to  which  the  damage  extends,  but  at  a  much  lower  rate  and 
with  no  noticeable  change  in  near-surface  yield.  The  anneals 
at  600, 700,  and  800  "C  for  30  s  can  be  seen  to  result  in  succes¬ 
sively  lower  near-surface  yields.  A  channeled  RBS  spectrum 
indistinguishable  from  the  unimplanted  spectrum  was 
achieved  after  800  *C  for  30  min  in  an  arsine  atmosphere 
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FIG  ?  Lattice  image  of  GaAs  following  amorphization  and  regrowth 
showing  microtwins  extending  to  the  surface 

Following  SPE  the  lowering  of  the  height  of  the  RBS  yield  at 
temperatures  of 200-600  ”C  is  believed  to  be  due  to  the  remo¬ 
val  of  dislocation  networks  and  of  faults  within  them.  A  var¬ 
iety  of  dislocation  network  sizes  and  configurations  results 
in  arrays  that  require  varying  temperatures  for  removal. 
This  would  account  for  a  near-surface  yield  that  is  tempera¬ 
ture  dependent.  Complete  removal  of  defects  accounting  for 
the  near-surface  yield  does  not  occur  unless  the  temperature 
is  increased  to  the  range  of  700  °C.  This  appears  to  be  due  to 
the  thermal  energy  required  for  removal  of  the  microtwins 
that  intersect  the  surface. 

In  conclusion,  a  correlation  of  BTE  calculations  with 
RBS  and  TEM  results  has  been  used  to  determine  the  frac¬ 
tional  displacement  of  Ga  and  As  atoms  necessary  for  the 
crystalline  to  amorphous  transition  in  GaAs  at  liquid-nitro¬ 
gen  temperature.  BTE  calculations  of  damage  profiles  are 
also  presented  and  correlated  with  stacking  faults  and  mi¬ 
crotwins  observed  in  TEM  results.  Finally,  RBS  data  for  a 
GaAs  anneal  sequence  are  correlated  with  TEM  observa¬ 
tions  to  identify  defects  responsible  for  the  RBS  yield  and  to 
determine  how  they  anneal.  Although  many  stacking  faults 
and  dislocations  are  removed  at  lower  temperatures,  tem¬ 
peratures  in  the  range  of  700  °C  are  required  for  complete 
microtwin  removal.  For  this  reason,  amorphization  fol¬ 
lowed  by  SPE  as  a  low-temperature  doping  technique  wi’l  be 
applicable  only  to  conditions  in  which  the  precipitation  of 
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dislocation  and  microtwin  nuclei  can  be  suppressed  as  was 
apparently  the  case  for  the  low-energy  shallow  implant  re¬ 
ported  previously. i: 
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ratory  program  through  the  Center  for  Materials  Research 
at  Stanford. 


'B.  L  Crowder.  J.  Electrochem.  Soc  117,  671  ( 1970) 

;J  S  Williams  and  M.  W  Austin.  Appl  Phys.  Lett.  36.  994  ( 1980) 

M  G.  Gnmaldi,  B.  M  Paine.  M.  A  Nicolet.  and  D  K.  Sadana.  J  Appl 
Phys.  52.  4032  (1981) 

JA  K  Rai,  R  S.  Bhattacharya.  P.  P  Pronko.  and  S.  C  Ling.  I  Appl  Phvs. 
54,2307  (1983). 

'D.  K  Sadana,  T.  Sands,  and  J.  Washburn.  Appl  Phvs.  Lett  44.  623 
(1984) 

Sands,  D.  K  Sadana.  R.  Gronsky,  and  J.  W'ashbum.  Appl.  Phys.  Lett. 
44,874  (1984). 

C.  Liccope,  Y.  I.  Nisstm,  and  C.  Meriadec,  J.  Appl.  Phys.  58.  3094  (1981). 
“L.  A  Chnstel  and  J  F.  Gibbons.  J  Appl  Phys.  52.  5050  ( 1981 ) 

T.  A  Christe),  J.  F.  Gibbons,  and  T  W  Sigmon.  J  Appl.  Phys  52.  7143 
(1981). 

'“M  D.  Giles  and  J  F  Gibbons,  IEEE  Trans  Electron  Devices  ED-32. 
1918  (1985). 

''J.  C  Bravman  and  R  Sinclair,  J  Electron  Microsc.  Tech.  1,  53  ( 1984) 
'•W.  G.  Opyd  and  J.  F.  Gibbons,  Mat  Res.  Soc  Symp  Proc  45,  273 
(1985). 


22 


Paper  #4 


Precipitation  of  Impurities  in  GaAs  Amorphized  by  Ion  Implantation 
W.G.  Opyd 

Lockheed  Research  and  Development  Division,  Palo  Alto,  California  94304-1 187 
J.F.  Gibbons 

Stanford  Electronics  Laboratories,  Stanford  University,  Stanford,  California  94305 
AJ.  Mardinly 

Lockheed  Missiles  and  Space  Company,  Sunnyvale,  California  94088-3504 

Abstract 

Impurities,  in  a  GaAs  layer  that  had  been  amorphized  by  ion  implantation, 
were  observed  to  precipitate  upon  annealing.  Photoluminescence  spectra  indicated 
that  the  resulting  high  electrical  resistivity  could  be  attributed  to  the  formation  of 
neutral  impurity  complexes  rather  than  a  compensation  mechanism.  Impurities 
studied  were  implanted  Si  and  Se.  Transmission  electron  microscopy  and  x-ray 
microanalysis  were  used  to  identify  impurity  precipitates  and  related  stacking  fault 
tetrahedra.  These  results  correlate  with  similar  examples  of  poor  activation  for 
impurities  in  GaAs  grown  by  low  temperature  molecular  beam  epitaxy. 
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There  has  been  interest  for  several  years  in  the  solid  phase  epitaxy  (SPE)  of 
a  GaAs  layer  amorphized  by  ion-implantation.  Interest  was  originally  stimulated  by 
the  discovery  that  SPE  at  a  relatively  low  temperature  (600°C  for  a  Si  crystal) 
permits  ion-implanted  dopants  to  be  incorporated  onto  substitutional  sites  in  a  Si 
crystal  [1].  For  the  compound  semiconductor  GaAs  a  low  temperature  process  for 
impurity  incorporation  is  attractive  because  it  provides  both  reduced  material 
decomposition  and  reduced  impurity  redistribution  during  annealing.  In  addition, 
an  impurity  implanted  into  a  preamorphized  layer  would  have  a  relatively  abrupt 
depth  profile,  without  the  deep  channeling  tail  typical  of  implants  into  crystals  [2]. 
SPE  of  GaAs  was  first  demonstrated  by  Williams  and  Austin  in  1980  [3]  and  has 
recently  been  investigated  by  others  [4-10].  Although  SPE  in  GaAs  has  been 
demonstrated  at  very  low  temperatures  [3]  (180°C),  extended  defects  are  often 
observed  in  the  implanted  layer  following  SPE.  The  extended  defects  have  been 
reported  to  anneal  at  higher  temperatures  [10]  but  no  significant  electrical 
activation  has  been  attained  [2,11],  In  this  letter,  the  results  of  an  investigation  into 
the  cause  of  poor  electrical  activation  are  presented,  indicating  that  amorphization 
by  ion  implantation  may  be  an  effective  means  of  producing  an  electrically  inactive 
layer  in  GaAs. 

The  approach  to  the  investigation  into  electrical  activation  was  to  implant 
known  n-type  impurities  into  a  GaAs  wafer,  half  of  which  had  been  preamorphized 
by  arsenic  implantation.  After  anneals  of  various  times,  temperatures  and  ambients. 
Hall  measurements  were  performed  for  electrical  activity,  secondary  ion  mass 
spectroscopy  (SIMS)  was  used  to  determine  chemical  depth  profiles, 
photoluminescence  (PL)  was  used  to  investigate  impurity  site  location,  transmission 
electron  microscopy  (TEM)  was  used  for  defect  structure  identification,  and  energy 
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dispersive  x-ray  microanalysis  (EDX)  was  used  for  analyzing  defect  structure 
chemical  composition. 

Arsenic  ions  were  implanted  into  undoped  semi-insulating  (SI)  liquid 
encapsulated  Czochralski  (LEC)  (100)  GaAs  to  a  dose  of  3xl014  cm'2  at  lOOkeV. 

m 

This  produced  an  amorphous  layer  of  800  A  as  confirmed  by  TEM.  Half  of  the 
wafer  was  masked  from  this  amorphizing  implant  but  not  from  the  subsequent 
impurity  implant  so  that  chemical  depth  profiles  and  electrical  activation  could  be 
compared  for  the  two  cases.  The  implants  were  intentionally  shallow  to  emphasize 
the  difference  in  the  two  halves  of  the  wafer  due  to  channeling  in  the  side  without 
the  arsenic  implant.  The  impurity  implants  tested  were  Si  at  25keV  and  Se  at 
50keV,  both  to  doses  of  1013  cm'2.  These  resulted  in  similar  chemical  depth  profiles 
with  SIMS  data  for  Si  shown  in  Fig.  1.  Rapid  thermal  anneals  (RTA)  were  used  for 
SPE  and  defect  removal  to  minimize  thermal  exposure.  Though  tests  have  been 
performed  with  arsenic  ambients  and  various  capping  layers,  the  data  presented  is 
primarily  for  an  ambient  of  forming  gas  with  a  silicon  proximity  cap. 

Electrical  activation,  as  measured  by  the  van  der  Pauw-Hall  technique, 
resulted  in  no  n-type  electrical  activity  for  the  Si  in  preamorphized  GaAs  following 
anneals  up  to  temperatures  at  which  conversion  to  p-type  is  expected  [12].  The 
percent  activation  as  a  function  of  anneal  temperature  for  a  RTA  of  5  s  is  shown  in 
Fig.  2  for  the  crystalline  and  amorphous  cases.  The  maximum  activation  for  the 
crystalline  case  is  low  as  expected  due  to  the  high  doping  level,  in  excess  of  1018  cm' 
To  see  if  the  lack  of  n-type  activity  was  due  to  compensating  defects  that  require 
anneal  temperatures  in  excess  of  the  temperature  at  which  type  conversion  occurs 
for  Si  doping,  Se  doping  was  attempted  since  it  is  known  to  activate  well  at  higher 
temperatures.  Though  type  conversion  was  not  observed  for  the  Se  implants,  the 
preamorphized  material  showed  no  activation  while  the  crystalline  control  activated 
as  expected  for  anneals  of  1000°C  to  1 100°C  for  3  s. 
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In  an  attempt  to  determine  the  site  location  of  the  implanted  Si,  PL 

o 

measurements  at  12°K  were  performed  for  wavelengths  of  8000  to  11000  A  with  a 
resolution  of  7  A.  The  PL  spectra  shown  in  Fig.  3  demonstrate  effects  in  the  region 
of  a  broad  luminescence  line  attributed  to  a  carbon  donor-acceptor  pair 
recombination  [13].  As  anneal  temperature  is  increased,  the  lines  for  the 
preamorphized  and  unimplanted  samples  show  no  significant  change,  while  the  Si 
implant  into  crystalline  GaAs  shows  a  line  developing  at  8350  A  that  has  been 
attributed  to  Si  on  an  As  site  [12,13].  This  represents  the  expected  shift  of  the  Si 
from  Ga  to  As  sites  which  has  been  used  to  explain  the  type  conversion  [12].  The 
absence  of  the  Si  impurity  lines,  however,  has  been  used  as  an  indication  that  the  Si 
is  complexed  into  an  electrically  inactive  defect  structure  and  not  simply 
compensated  by  impurities  of  the  opposite  type[13]. 

In  order  to  identify  the  defect  structure  that  may  be  preventing  the  impurities 
implanted  into  preamorphized  GaAs  from  activating,  high  resolution  TEM  was 
performed  on  cross-section  specimens.  A  TEM  micrograph  of  the  GaAs  viewed  in  a 
<  1 10  >  direction  and  with  lattice  resolution  is  shown  in  Fig.  4  for  a  specimen 
implanted  with  Si  and  annealed  at  1050°C  for  3  s.  Defects  identified  were  stacking 
fault  tetrahedra  (SFT)  and  precipitates.  Under  the  diffracting  conditions  of  Fig.  4, 
these  two  structures  may  appear  distinct,  but  upon  closer  analysis  each  precipitate 
was  found  to  be  associated  with  a  SFT.  The  pricipitates  were  confirmed  to  contain 
Si  by  EDX  measurements.  Identical  structures  were  noted  for  specimens  in  which 
Se  had  been  implanted  into  the  preamorphized  GaAs.  Precipitates  at  the  core  of 
SFTs  have  previously  been  noted  to  correlate  with  poor  electrical  activation  and 
were  attributed  to  local  stoiciometric  disturbances  caused  by  an  implant  [14]. 

These  results,  indicating  a  lack  of  Si  activation,  correlate  weil  with 
experiments  on  low  temperature  molecular  beam  epitaxy  (LT-MBE)  [15].  In  .he 
LT-MBE  case,  activation  of  Si  incorporated  during  growth  was  assumed  dependent 


upon  the  growth  rate  being  slow  enough  to  allow  sufficient  surface  migration  of  the 
Ga  and  As  to  prevent  local  stoichiometric  disturbances.  If  the  growth  rate  was  too 
fast,  lack  of  activation  was  assumed  to  be  due  to  compensation  by  acceptors  but  no 
evidence  for  this  was  given.  A  comparison  of  the  rate  for  SPE  [8]  with  the  maximum 
LT-MBE  rate  for  Si  activation  [15]  is  shown  in  Fig.  5.  Though  this  comparison  can 
not  be  made  directly,  since  the  LT-MBE  is  limited  by  surface  migration  while  SPE  is 
limited  by  bulk  diffusion,  it  implies  that  activation  of  impurities  in  GaAs  following 
SPE  might  not  be  expected. 

In  cor  'lusion,  impurities  implanted  into  preamorphized  GaAs  have  been 
shown  to  precipitate,  inhibiting  electrical  activation.  Though  this  conclusion  rules 
out  the  use  of  implants  into  preamorphized  GaAs  for  achieving  abrupt  doping 
profiles,  it  suggests  a  technique  for  producing  electrically  inactive  material  in  spite 
of  the  presence  of  impurities.  Amorphization  of  a  layer  of  GaAs  containing 
impurities  would  cause  the  local  stoichiometric  disturbances  necessary  to  produce 
the  defects  that  collapse  into  SFTs  and  precipitate  the  impurities. 

We  would  like  to  thank  LJ.  Dries  and  S.K.  Ichiki  of  Lockheed  for  assistance 
in  PL  and  TEM  work  and  W.S.  Lee  and  J.S.  Harris  of  Stanford  for  assistance  in 
correlation  with  LT-MBE.  This  work  was  supported  at  Lockheed  by  independent 
research  funds  and  at  Stanford  by  the  Army  Research  Office  under  contract 
DAAG29-85-K-0054. 
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Figure  Captions 


1.  Depth  profile  of  implanted  Si  as  determined  by  SIMS  for  implants  into  crystalline 
and  amorphous  GaAs. 

2.  Percent  electrical  activation  of  101-3  Si/cm2  implanted  at  25keV  into  crystalline 
and  amorphous  GaAs  as  a  function  of  temperature  for  5  s  anneals. 

3.  Photoluminescence  spectra  of  GaAs  demonstrating  the  development  of  the  Si  on 
As  site  line  with  increasing  anneal  temperature  for  implants  into  crystalline  material 
while  this  line  is  not  present  in  preamorphized  or  unimplanted  material. 

4.  High  resolution  TEM  micrograph  of  Si  implanted,  preamorphized  GaAs  after  an 
anneal  of  1050°C  for  3  s. 

5.  Comparison  of  the  rate  for  SPE  with  the  maximum  rate  of  LT-MBE  for  which  Si 
activation  has  been  observed. 
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We  have  demonstrated  a  new  technique  for  III-V  epitaxial  layer  growth  combining  rapid 
thermal  processing  and  metalorganic  chemical  vapor  deposition.  This  technique  yields 
enhanced  layer  thickness  control  and  abrupt  interfaces  while  maintaining  a  high  growth  rate 
(>  10  A/s).  Multilayer  structures  have  been  grown  with  smooth,  featureless  surfaces  and  good 
electrical  quality  (Nd  =  2X  10l6cm~\/i„  =  3000  cmW  s)  using  trimethylarsenic  and 
trimethylgallium. 


The  recent  interest  in  GaAs  and  related  III-V  materials 
has  been  stimulated  by  the  development  of  epitaxial  growth 
techniques  which  yield  materials  of  both  good  crystal  and 
electrical  quality  along  with  tight  control  over  layer  thick¬ 
ness  and  interface  abruptness.  One  of  these  techniques,  me¬ 
talorganic  chemical  vapor  deposition  ( MOCVD ) ,  produces 
abrupt  interfaces  using  elaborate  and  costly  gas  switching 
systems  along  with  reduced  reactor  pressure.  Recently,  a 
new  method  of  layer  control  has  been  demonstrated  by  Gib¬ 
bons  et  al.  for  high  quality  thin  layer  growth  of  silicon  layers 
with  abrupt  interfaces. ' 2  This  method,  limited  reaction  pro¬ 
cessing  (LRP),  uses  rapid  precise  changes  in  the  substrate 
temperature  to  switch  layer  growth  on  and  off  rather  than 
relying  on  gas  phase  switching. 

We  report  rapid  thermal  switched  growth  of  GaAs  epi¬ 
taxial  layers.  We  have  extended  this  technique  with  equal 
success  to  AI^Ga,  _,As  (0<x<0.25)  and  lnvGa,_,As 
( 0  <  y  <  0.03 ) ,  all  on  GaAs  substrates.  Good  electrical  char¬ 
acteristics  and  surface  morphology,  featureless  to  within  the 
resolution  of  our  Nomarski  microscope  ( 1 100  X  ),have  been 
consistently  obtained.  In  contrast  to  standard  MOCVD.  this 
technique  maintains  a  high  growth  rate  ( 10  A/s)  at  atmo¬ 
spheric  pressure  while  producing  abrupt  interfaces. 

Layers  were  grown  using  trimethylarsenic,  trimethyl¬ 
gallium,  trimethylaluminum,  and  tnethylindium,  all  pur¬ 
chased  from  Alfa  Products  (Danvers,  M A).  Substrates,  un¬ 
doped  and  Si  doped  ( 100)  Czochralski,  were  degreased  and 
given  a  5:1:1  H:S04:H:0::H:0  etch  prior  to  loading.  Opti¬ 
mum  layers  were  grown  at  670-720  °C  with  3.5  l/min  of 


total  H,  flow  and  a  trimethylarsenic  partial  pressure  of  about 
0.5  Torr.  GaAs  layers  have  n-type  background  doping 
(2Xl016  cm“3<Arrf<lXl017  cm-3)  and  mobilities  of 
2500-3000  cm2/V  s  at  room  temperature.  We  believe  the 
high  backgrounds  are  due  to  impurity  incorporation  from 
the  trimethylarsenic,  as  has  been  reported  by  others.3  Tri¬ 
methylarsenic  was  selected  on  the  basis  of  safety;  it  is  known 
that  arsine  produces  much  higher  purity  films. 

Our  reactor  design  is  shown  in  Fig.  1 .  Since  layer  switch¬ 
ing  is  now  done  thermally,  it  is  possible  to  use  a  simplified 
gas  control  system.  The  wafer  sits  on  a  thin  ( 20  mil )  graphite 
susceptor  heated  from  the  underside  by  a  bank  of  high  power 
tungsten  halogen  lamps,  as  in  a  rapid  thermal  annealer. 
Temperature  is  measured  with  a  thermocouple  inserted 
down  a  2-mm-o.d.  sealed  quartz  tube  in  contact  with  the 
susceptor.  For  initial  calibration  an  additional  thermocou¬ 
ple  is  welded  to  the  sample4  and  the  relationship  between  the 
wafer  temperature  and  the  inserted  thermocouple  is  estab¬ 
lished.  The  graphite  susceptor  is  used  for  temperature  mea¬ 
surement  only,  we  believe  that  supporting  the  sample  on 
quartz  pins  alone  and  measuring  temperature  with  an  opti¬ 
cal  pyrometer  would  improve  this  technique  by  reducing 
thermal  rise  and  fall  times. 

In  rapid  thermal  MOCVD  gas  flows  are  initiated  and 
stabilized  while  the  wafer  is  cool.  The  growth  of  epitaxial 
layers  is  initiated  by  pulsing  the  lamps  and  bringing  the  wa¬ 
fer  rapidly  to  growth  temperature.  At  the  end  of  each  layer 
the  lamps  are  shut  off,  the  wafer  cools  rapidly  towaru  room 
temperature,  and  the  reaction  is  halted.  The  cooling  rate  is 
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FIG.  I .  Schematic  drawing  of  the  quartz  reac¬ 
tion  chamber  used  in  our  experiments  The 
wafer  sits  on  a  low  thermal  mass  graphite  sus¬ 
ceptor  and  temperature  is  monitored  by  a 
thermocouple  sheathed  in  a  2-mm-o.d  quartz 
tube. 
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enhanced  by  cooling  the  walls  of  the  quartz  reaction 
chamber  with  a  high  flow  of  compressed  air.  switched  on  at 
the  termination  of  layer  growth.  Heat  conducts  rapidly  away 
from  the  wafer  through  the  hydrogen  carrier  gas.  Rise  and 
fall  times  are  on  the  order  of  10  s.  Any  desired  structure  may 
be  grown  by  a  sequence  of  these  steps,  as  shown  in  Fig.  2  fora 
GaAs/AI,  Ga,  _  ,  As  multilayer  structure.  Note  that  the 
GaAs  is  grown  at  670  °C  and  the  AiGaAs  at  720  °C.  We 
found  these  temperatures  gave  the  best  mobility  and  surface 
morphology  for  single  layers  of  each  material. 

One  of  our  early  concerns  w  as  that  poor  material  w  ould 
be  grown  during  the  cool  down  period,  creating  defects  in 
the  growth  of  subsequent  layers.  Test  structures  consisting 
of  multiple  layers  of  GaAs  or  alternating  layers  of  GaAs/ 
AiGaAs  (3-11  layers.  500-2500  A  layer  thickness)  have 
been  grown  by  a  sequence  of  rapid  thermal  cycles.  These 
structures  have  excellent  surface  morphology  and  electrical 
characteristics,  equivalent  to  single  layers  grown  by  conven¬ 
tional  MOCVD  in  our  reactor.  Such  test  structures  have 
been  analyzed  by  Rutherford  backscattering  and  ion  chan¬ 
neling,  showing  no  evidence  of  crystal  defects  at  ihe  inter¬ 
faces.  Minimum  yields  are  about  4.07r,  the  same  as  bare 
GaAs  wafers.  Rutherford  backscattering  detects  only  gross 
crystal  defects;  more  sensitive  tests  of  interface  quality  are  in 
progress,  including  device  fabrication  and  high  resolution 
transmission  electron  microscopy. 

An  Auger  depth  profile  of  a  GaAs/Al,,  ,.Ga„  ..As  mul¬ 
tilayer  structure  grown  by  this  technique  is  shown  in  Fig.  J. 
Layers  are  abrupt  to  within  the  resolution  of  Auger  (  50-00 
A).  The  total  thickness  of  this  structure  is  about  fiOOO  A. 
Because  gases  are  completely  purged  between  layers,  we  be¬ 
lieve  this  technique  is  capable  of  atomically  abrupt  inter¬ 
faces. 

In  conclusion,  we  have  demonstrated  mpid  thermal 
MOCVD  hy  growing  multiple  layers  of  lucn  uualily  III-V 


compound  semiconductors  These  layers  exhibit  electrical 
characteristics  and  surface  morphology  equivalent  to  layers 
produced  by  traditional  MOCVD  in  our  reactor.  The  rapid 
thermal  MOCVD  technique  promises  abrupt  interfaces  and 
thin  layers  combined  with  high  growth  rates. 

W’e  would  like  to  acknowledge  the  technical  assistance 
of  K.  E.  Williams  and  helpful  discussions  with  D  Houng  of 
Hewlett-Packard  Laboratories.  C.  P  Kuoof  Hewlett-Pack¬ 
ard  Optoelectronics  Division,  and  H  Beneking  of  Aachen 


•35  3EPth  profile 

100  - — - - 


80  - 


tU 

O  60  •  . 

C  '  ,  «  , 

OJ  V.'v- 

O.  -  ,  ■  •  ,  > 

U  1  NA,  'y-V*  v*v- 


>PL;T,rEP«?4G  ’  VE 


M(i  •  SpulterinL' Aucer  pr*'fiie  Tor  ".imr't*  »  ,«howin  !  a  v.fruv  •  '  •  ■' 

,\'ri!.iiii;k’  iuA\  \l  li.i _  \s  ,.i\crv  T  :al  Thickness  i'  jdoi \ 

Xiurmtmm  scjlcil  h\  2  (  hjfii's  loans  jnJ  i.iics.  KcJw  1  » 

(.  \  i 


35 


Paper  #5 


Technical  Institute.  Aachen.  West  Germain  We  also  ac¬ 
knowledge  the  support  ot'J  Z.nuda  and  D  A  Reynolds  at 
Defense  Ads anced  Research  Protects  Agency  (contracts 
DAAG2‘)-S5-K-(KJ54  and  DAAG24-k5-K  2>'>.  S  Reyn¬ 
olds  acknowledges  the  support  of  an  Office  of  Nasal  Re¬ 
search  fellowship. 


J  I  Gibbons. (_  M  ( iioiicl  .and  K  f  ^  lihaips.  A ppi  P!i>s  l.cH  4""  '2! 
;  t‘^5) 

C  VI  Gronei.  J  C  Siurm.  K.  E  G  illums,  <  F  Gibbons,  and  S  D  A  *J- 
soii.  AppE  I’li\\  Ecu  48.  I  * 1 1 2  !^>n- 

C  A  IroniMiu.  I1  Gih.iri.  K  Dnnihe.  \  MoniviE  J  IJouu.  and  h  EJ 
Jam,  Re\  Phys  Appl  i  France)  20.  fb'J  <  !4*5) 

L  Hoyf  and  J  F  Gibbons.  Maier  Res  St*  Ssmp  Pr<>c  52.  15  i  N“bj 


Rapid  thermal  annealing  of  Si-implanted  GaAs  with  trimethylarsenic 
overpressure 
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We  have  developed  a  novel  rapid  thermal  processor  to  perform  annealing  of  ion-implanted 
GaAs  in  a  trimethylarsenic  overpressure.  This  has  allowed  study  of  arsenic  ambient  annealing 
in  a  time/temperature  regime  not  accessible  with  an  arsine  furnace.  We  have  compared  Si 
implant  activation  efficiency  and  surface  degradation  for  arsenic  ambient  and  proximity 
capped  anneals.  The  arsenic  ambient  gives  consistently  higher  activation  efficiency  with  better 
surfaces. 


The  annealing  of  ion-implanted  GaAs  has  received 
much  study  because  of  its  practical  importance  for  field-ef¬ 
fect  transistor  fabrication.1  Preferential  As  evaporation 
from  the  surface  must  be  avoided.  This  is  normally  done 
either  by  providing  an  arsenic  overpressure  in  a  hot  wall 
furnace,  or  by  rapid  thermal  annealing  (RTA)  with  a  depos¬ 
ited  encapsulated  or  proximity  cap.  For  practical  reasons, 
none  of  these  methods  is  entirely  satisfactory.  Arsenic  fur¬ 
naces  involve  long  thermal  exposure  and  use  large  quantities 
of  arsine.  Deposited  encapsulants  (such  as  Si,N4  or  SiO:) 
stress  the  surface  and  tend  to  crack,  causing  localized  in¬ 
creases  in  surface  deterioration. : 

We  report  rapid  thermal  annealing  of  GaAs  with  tri¬ 
methylarsenic  (TMAs)  overpressure,  which  has  advantages 
over  the  previous  techniques.  Thermal  exposure  can  be  mini¬ 
mized  and  arsenic  overpressure  need  only  be  supplied  for  a 
short  time.  Surface  morphology  is  better  preserved  when  no 
encapsulant  or  proximity  cap  touches  the  GaAs  surface.  As 
ambient  anneals  have  been  studied  for  short  times  unobtain¬ 
able  in  an  arsine  furnace,  and  As  ambient  RTA  has  been 
compared  with  the  proximity  capping  method.  We  have  also 
done  preliminary  experiments  with  nitride  encapsulation. 
Proximity  capping  was  used  for  comparison  because  it  is  a 
readily  duplicated  method.  Results  obtained  with  nitride 
caps  are  more  dependent  on  the  deposition  conditions,  so 
comparisons  between  laboratories  are  more  difficult. 

A  schematic  drawing  of  our  anneaier  is  shown  in  Fig.  1 . 
The  design  is  chosen  to  meet  the  problem  of  maintaining  an 
arsenic  overpressure  in  a  cold  wall  reactor.  A  very  thin  ( 2 
mil)  molybdenum  sheet,  rolled  into  a  semicylinder,  slides 
inside  a  quartz  tube  and  is  open  at  the  bottom  toward  the 
tungsten  halogen  lamps.  It  absorbs  enough  heat  to  prevent 
arsenic  depletion  near  the  wafer,  reaching  over  500  °C  about 
5  s  after  the  lamps  are  switched  on.  Without  it.  arsenic  con¬ 
denses  on  the  quartz  wails  and  reduces  the  effective  arsenic 
overpressure.  The  molybdenum  sheet  is.  in  effect,  a  rapid 
thermal  wall.  The  samples  are  supported  and  heated  indi¬ 
rectly  by  a  thin  (JO  mil)  graphite  susceptor.  A  1/64-in. 
stainless-steel  sheathed  thermocouple  is  inserted  in  the  sus¬ 
ceptor  for  closed  loop  temperature  control.  The  tube  is 
cooled  with  a  high  How  of  compressed  air,  switched  on  at  the 
end  ofthcanneji  cycle  Rive  and  fall  times  are  on  the  order  of 
10  s. 


Anneals  were  performed  in  1  1pm  of  flowing  H-  gas. 
with  or  without  TMAs.  TMAs  was  transported  in  vapor 
form  by  passing  H:  through  a  bubbler  held  at  —  5  *C.  A 
TMAs  partial  pressure  of  2  Torr  was  found  to  be  adequate  to 
protect  the  surface  up  to  850  °C  For  experiments  involving 
higher  temperatures,  a  pressure  of  4  Torr  was  used.  TMAs 
was  selected  because  it  was  less  toxic  and  easier  to  handle 
than  arsine.  In  this  application  the  additional  methyl  groups 
from  TMAs  do  not  appear  to  be  much  of  a  disadvantage.  Our 
rapid  thermal  annealer  is  cleaner  than  an  arsine  furnace  be¬ 
cause  TMAs  only  decomposes  while  the  susceptor  is  hot. 
Since  anneal  times  are  so  short,  comparatively  little  arsenic 
dust  is  generated. 

We  first  investigated  the  effect  of  annealing  time  on  im¬ 
plant  activation  efficiency.  We  made  Hall  effect  measure¬ 
ments  in  the  van  der  Pauw  configuration  to  determine  the 
sheet  earner  concentration,  then  divided  this  number  by  the 
implant  dose.  For  temperatures  j  8 JO  =C  anneal  time  has  lit¬ 
tle  effect  on  activation,  as  show  n  in  Fig.  2  for  a  10'-'  cm  \  50 
keV  Si  implant.  A  10-s  anneai  is  just  as  effective  as  a  much 
longer  one  for  implant  activation.  This  suggests  that  most  of 
the  time  a  wafer  spends  in  an  arsine  furnace  may  be  unneces¬ 
sary  thermal  exposure. 

Activation  efficiency  as  a  function  of  temperature  has 
been  studied  for  a  variety  of  Si  impiant  doses,  comparing 
TMAs  ambients  to  Si  proximity  caps.  Doses  ranged  from 
10"  to  10"  cm  -  In  ail  cases  the  highest  activation  effi¬ 
ciency  was  obtained  with  TMAs  overpressure.  Resuits  for  a 
101 1  cm  ~  :,  50  keV  implant  are  shown  in  Fig.  J.  Optimum 


Fl(i  1  Sc  hem  j  tic  draw  me  of  ‘he  rapid  --.ernai  annealer  u-vrd  for  :ri* 
neih\ lursenic  o\ erpressure  annexing. 


37 


mi;ian*  S 


Time  (seconds) 

FIG  2  Activation  efficiency  vs  anneal  time  for  TM  As  ambient  anneals 

annealing  temperatures  were  870  °C  for  the  proximity  cap 
and  900  °C  for  TMAs,  yielding  electron  mobilities  of  2600 
cnr/V  s  in  both  cases.  The  corresponding  activations  were 
50-55%  for  the  proximity  cap  and  60-65%  for  the  TMAs. 
The  improvement  afforded  by  TMAs  increases  at  higher  an¬ 
nealing  terroeratures.  Secondary  ion  mass  spectroscopy 
(SIMS)  denth  profiling  verified  that  no  Si  diffusion  oc¬ 
curred  for  10  s  anneals  at  900  °C.  in  TMAs  or  with  Si  proxim¬ 
ity  caps. 

Some  of  our  implant  samples  were  encapsulated  with 
Si,N,  by  a  low-pressure  plasma  enhanced  process  at  200  °C. 
The  maximum  activation  efficiencies  were  about  two-thirds 
as  high  as  obtained  for  similar  samples  annealed  in  TMAs. 
Visible  blistering  of  the  nitride  was  observed  for  10  s  anneals 
ut  900  °C  and  above.  After  removing  the  nitride,  damage  to 
the  GaAs  surface  of  these  samples  was  observed  under  our 
Nomarski  microscope  ( unannealed  samples  with  the  nitride 
removed  had  no  damage ) .  Surface  damage  was  occasionally 
observed  even  for  those  lower  temperature  samples  which 
did  not  blister.  Other  workers  have  obtained  good  results 
with  nitride  encapsulation,  so  our  results  with  it  should  not 
be  taken  as  conclusive. Nitride  encapsulation  seems  to  be 
a  process  which  must  be  optimized  by  each  laboratory  using 
it. 
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sure  .done  remained  featureless  lor  ,.il  .umeui  -  rijn:.  j;> 
file  surface  protected  h\  a  Si  proximity  cap  alone  wilho  I 
TMAs  began  showing  signs  of  localized  pitting  h  r  10  s  an¬ 
neals  at  925  °C.  The  surface  degradation  was  visible  by  eve 
for  10  s  anneals  at  950 *C.  The  silicon  proximny  caps  were 
inspected  after  each  anneal  and  replaced  it'  i  here  was  any  sign 
of  discoloration. 

In  our  first  experiments.  GaAs  waters  were  used  as 
proximity  caps.  We  unexpectedly  discovered  that  GaAs 
proximity  caps  are  much  less  effective  than  Si  caps  at  pro¬ 
tecting  the  surface.  Despite  the  use  of  a  new  GaAs  cap  for 
each  anneal,  the  proximity  caps  themselves  seem  to  decom¬ 
pose  and  promote  further  decomposition  of  the  underlying 
wafer.  With  GaAs  caps,  surface  decomposition  v  isible  lo  the 
eye  begins  as  low  as  850  °C  for  a  10-s  anneal. 

TMAs  overpressure  is  thus  much  more  effective  than 
proximity  capping  at  preventing  gross  surface  decomposi¬ 
tion.  Auger  sputter  profiling  was  used  to  make  a  qualitative 
comparison  of  surface  composition  at  lower  temperatures, 
where  both  As  overpressure  and  proximity  capping  yielded 
featureless  surfaces  under  the  microscope.  The  differences 
revealed  by  Auger  were  smaller  than  expected  (for  830  °C 
anneals).  Figured  shows  a  set  of  Auger  profiles  for  a  sample 
annealed  in  TMAs  ( 10  s  at  830  °C).  then  exposed  to  room 
air.  The  most  notable  feature  is  the  oxy  gen  signal  extending 
about  120  A  into  the  sample.  A  decrease  in  the  As  signal 
extends  to  about  the  same  depth.  Samples  annealed  with  a 
proximity  cap  are  very  similar,  bui  appear  to  develop  a 
slightly  thicker  native  oxide.  To  quantitatively  compare  the 
extent  of  oxidation,  we  integrated  (he  oxygen  and  arsenic 
signals  lo  determine  the  relative  O  gain  and  As  loss.  The 
results  are  shown  in  Table  I  While  noth  anneai  conditions 
generate  some  surface  deterioration.  TMAs  annealed  sur¬ 
faces  are  better. 

In  summary,  we  have  designed  a  TMAs  ambient  rapid 
thermal  annealer  which  maintains  an  As  overpressure  in  a 
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I  ABLE  I  Comparison  of  ilie  o\\  gen  pain  and  An  loss  lor  samples  annealed 
under  various  conditions  (and  (hen  exposed  loumm  air),  obtained  tn  mte- 
gratme  the  Auger  oxvcen  and  arsenic  signals 


Sample 

Annealing 

conditions 

O  gained  tn 
surface  oxide 
( 101'  atoms/cm 

As  IcHi  from  bulk  lo 
ambient  or  into 
surface  oxide 
( I01'  aioms/crrr ) 

i 

as-implanted 

unannealed 

5.1 

2.5 

34QI 

10  s  @ 

830  ‘ C 
TMAs 

12.1 

rs 

3431 

10  s  @ 
830*C 

Si  prox  cap 

15.7 

8.7 

nominally  cold  wall  reactor.  These  are  the  first  RTA  experi¬ 
ments  using  a  controlled  As  ov  erpressure  similar  to  an  arsine 
furnace.  We  have  observed  that  short  anneals  in  As  appear  to 


be  as  effective  as  long  (Mies  for  activating  implants,  -o  that 
ihc  added  thermal  exposure  of  an  arsine  furnace  nuv  be  un¬ 
necessary.  Rapid  As  ambient  anneals  give  higher  activation 
of  Si  implants  and  less  surface  decomposition  than  proximity 
capping,  as  expected. 
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acknowledges  the  support  of  an  Office  of  Naval  Research 
Fellowship.  We  thank  S.  Bui  of  Avantek  for  technical  assis¬ 
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Open-tube  Zn  diffusion  in  GaAs  using  diethylzinc  and  trimethylarsenic: 
Experiment  and  model 
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We  have  characterized  the  diffusion  of  Zn  into  GaAs  from  the  organometallic  sources 
diethylzinc  and  trimethylarsenic.  This  method  produces  surface  hole  concentrations  in  excess 
of  10io  cm  -  3  with  good  control  of  junction  depths  as  shallow  as  0. 1  pm.  Smooth  surface 
morphology  is  retained.  The  profile  shape  is  much  more  complex  than  the  accepted  interstitial- 
substitutional  Zn-diifusion  model  would  predict.  To  explain  the  observed  profiles,  a  new  model 
for  Zn  diffusion  is  proposed  and  implemented  in  a  computer  simulation. 


I.  INTRODUCTION 

Zn  diffusion  in  GaAs  is  an  important  processing  step  in 
the  fabrication  of  diodes,  junction  field-effect  transistors 
(JFETs),  and  bipolar  transistors.1  Unfortunately,  control 
over  the  diffusion  process  has  often  been  difficult  and  can 
limit  device  performance.  Widely  varying  results  with  Zn 
diffusion  have  been  reported.  An  understanding  of  the  phys¬ 
ical  mechanisms  of  the  diffusion  process  may  allow  it  to  be 
brought  under  better  control. 

Vapor  phase  diffusions  using  a  diethylzinc  (DEZn) 
source  were  reported  in  a  paper  by  Dohsen  et  air  These 
authors  fabricated  GaAs  JFETs  by  localized  diffusion 
through. a  patterned  nitride  (Si3N4)  mask.  We  employed 
this  vapor  phase  diffusion  method  to  perform  blanket  diffu¬ 
sions  for  the  fabrication  of  p* -i-n  diodes.  In  doing  so.  we 
observed  that  the  junction  depth  and  profile  shape  are  a  com¬ 
plicated  function  of  the  diffusion  conditions.  We  report  our 
experimental  results  using  this  method,  and  we  propose  a 
new  model  to  explain  the  complex  profile  shapes. 

It  is  widely  accepted  that  Zn  diffuses  in  GaAs  by  an 
interstitial-substitutional  mechanism,  similar  to  that  pro¬ 
posed  by  Frank  and  Turnbull3  for  Cu  in  Ge.  In  this  model,  a 
small  fraction  of  the  total  dissolved  Zn  is  present  interstitial- 
ly.  This  interstitial  Zn  diffuses  rapidly  and  dominates  the 
diffusion  process.  This  model  predicts  an  abrupt  concentra¬ 
tion-dependent  diffusion  profile  where  the  effective  diffusion 
coefficient  is  proportional  to  the  square  of  the  total  Zn  con¬ 
centration.4  This  model  works  well  in  the  high-concentra¬ 
tion  surface  region  of  our  Zn  profiles,  but  it  does  not  predict 
the  complex  double  profiles  that  we  and  others5  have  ob¬ 
served.  Our  model  is  a  modification  of  the  interstitial-substi¬ 
tutional  mechanism. 

II.  EXPERIMENTAL  PROCEDURE 

Diffusions  were  carried  out  in  a  hot-wall  quartz  reaction 
chamber  (40  mm  diam)  at  atmospheric  pressure,  with  1  1/ 
min  of  H2  earner  gas  flow.  Trimethylarsenic  (TMAs)  and 
DEZn  were  transported  in  vapor  form  by  flowing  H2 
through  bubblers  held  at  constant  temperature.  The  TMAs 
and  DEZn  partial  pressures  (pTMA,  and  pDlz„ )  were  varied 
over  a  range  from  0. 1  to  3.0  Torr.  TMAs  was  used  to  prevent 
GaAs  surface  decomposition  at  elevated  temperatures. 

GaAs  wafers  were  undoped  ( 100)  Czochralski  material 


polished  on  one  side,  obtained  from  Morgan  semiconductor. 
The  dislocation  densities,  as  specified  by  the  manufacturer, 
ranged  from  2x  104  to  4.3 X  !04  cm-2.  Samples  were  de¬ 
greased  with  organic  solvents  prior  to  loading.  When  we  first 
observed  the  complex  profile  shape,  we  wondered  if  it  was 
due  to  residua]  surface  damage  from  the  manufacturer.  One 
set  of  samples  was  degreased  and  etched  for  5  min  in  5:1:1 
H2S04:H202:Hj0  prior  to  loading.  By  simultaneously  pro¬ 
cessing  another  set  of  samples  which  were  degreased  only, 
we  verified  that  etching  the  wafers  had  no  effect  on  diffusion 
rates  or  profile  shape.  To  be  consistent,  all  samples  included 
in  the  data  presented  here  were  degreased  only. 

Samples  were  loaded  into  the  reaction  chamber  on  a 
graphite  susceptor.  The  chamber  was  evacuated,  then  back¬ 
filled  with  H2  and  purged.  TMAs  flow  was  started  and  the 
sample  was  ramped  to  the  diffusion  temperature  in  3  min. 
After  1  min  at  the  diffusion  temperature  the  DEZn  was 
switched  on.  At  the  end  of  the  diffusion  time  the  DEZn  was 
switched  off,  the  chamber  was  purged  for  1  min.  then  cooled. 
When  the  susceptor  reached  500  °C  the  TMAs  was  switched 
off.  Total  cooling  time  to  below  100  °C  was  less  than  5  min. 

Controlling  the  diffusion  by  switching  the  DEZn  flow 
was  found  to  be  the  most  satisfactory  method.  If  the  suscep¬ 
tor  is  cooled  rapidly  with  DEZn  flowing,  elemental  Zn  con¬ 
denses  on  the  wafer.  Thus,  temperature  switching  is  ineffec¬ 
tive.  By  switching  on  the  DEZn  flow  for  a  period  of  several 
tens  of  seconds,  diffusions  could  be  well  controlled  to  depths 
as  shallow  as  0.1  pm.  Well-controlled  shallower  diffusions 
might  be  obtained  by  increasing  the  H2  flow,  thus  purging 
the  chamber  more  quickly.  The  only  disadvantage  of  switch¬ 
ing  the  DEZn  flow  to  control  junction  depth  is  a  small 
amount  of  out-diffusion  near  the  GaAs  surface,  which  can  be 
seen  in  Fig.  1.  Tuck  and  Houghton6  showed  that  when  Zn- 
doped  GaAs  is  annealed  without  a  Zn  source,  very  little 
further  movement  of  the  Zn  profile  occurs.  Thus.  Zn  diffuses 
rapidly  only  when  it  is  being  incorporated  at  the  surface.  Our 
results,  discussed  in  Sec.  Ill  confirm  this  observation.  This 
behavior  is  important  in  developing  a  Zn-difTusion  model. 

Samples  diffused  over  the  temperature  range  from  650- 
750  'C  had  smooth,  specular,  and  generally  featureless  sur¬ 
faces  when  examined  under  our  Normarski  microscope. 
Samples  diffused  much  above  750  'C  begin  to  show  signs  of 
surface  roughening.  Samples  diffused  below  625  "C  show 
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FIG.  I.  Electrical  and  SIMS  profiles  for  a  5-min  diffusion  at  725  *C.  The 
surface  peak  is  a  SIMS  artifact,  but  the  dip  immediately  behind  it  is  prob¬ 
ably  a  small  amount  of  out-diffusion  from  the  1-min  final  purge  (see  discus¬ 
sion  in  Sec.  II). 


areas  of  what  we  believe  to  be  metallic  Zn  condensation  if 
high  DEZn  flows  are  used. 

The  walls  of  the  reaction  chamber  receive  fairly  heavy 
deposition  of  As,  Zn,  and  Zn,  As.  compounds,  comparable 
in  degree  to  a  GaAs  metalorganic  chemical  vapor  deposition 
(MOCVD)  reactor.  Since  this  is  a  hot-wall  reactor,  these 
deposits  can  have  an  unpredictable  effect  on  subsequent  dif¬ 
fusions  at  a  different  temperature  or  DEZn  flow.  Therefore, 
the  reaction  chamber  was  baked  out  at  a  high  temperature 
(  >900°C)  with  flowing  TMAs  between  each  diffusion  to 
remove  volatile  Zn  compounds  (it  was  removed  and  chemi¬ 
cally  cleaned  every  3-4  runs).  With  the  reaction  chamber 
baked  out  or  chemically  cleaned  between  runs,  we  were  able 
to  reproduce  junction  depths  to  w  ithin  =  15%.  We  verified 
that  samples  annealed  with  no  DEZn  flow  (TMAs  only) 
after  the  chamber  had  been  baked  out  showed  no  measurea- 
ble  conductivity.  Note  that  it  would  be  difficult  to  carry  out 
this  diffusion  method  in  a  cold-wall  reactor.  Zn  and  As 
would  diffuse  from  the  hot  zone  near  the  susceptor  and  depo¬ 
sit  on  the  cold  walls,  making  it  difficult  to  maintain  the  re¬ 
quired  Zn  and  As  partial  pressures. 

Zn  concentration  profiles  in  the  diffused  samples  were 
measured  by  secondary  ion  mass  spectrometry  (SIMS).7 
Electrical  measurements  were  made  by  the  differential  van 
der  Pauw  technique.  In  this  method,  sheet  resistance  and 
Hall  voltage  are  measured  at  the  GaAs  surface.  Carrier  con¬ 
centrations  are  obtained  by  successively  etching  away  thin 
GaAs  layers  and  then  remeasuring  sheet  resistance  and  Hall 
voltage.  We  have  analyzed  more  than  40  profiles  by  either 
SIMS,  electrical  measurements,  or  both 


III.  EXPERIMENTAL  RESULTS  lJaper  ^  7 

Figure  1  shows  electrical  and  chemical  profiles  for  a  5- 
min  diffusion  at  725  *C.  The  most  notable  feature  is  the  lone 
tail  extending  beyond  the  high-concentration  surface  region. 
Electrical  measurements  imply  that  virtually  all  of  the  Zn  is 
electrically  active,  probably  even  in  the  tail  region.  Electncal 
measurements  do  not  extend  to  the  same  depth  as  SIMS 
measurements  because  the  contacts  on  a  van  der  Pauw  sam¬ 
ple  generally  fail  before  the  conductive  layer  is  completely 
removed. 

Figure  2  shows  SIMS  profiles  for  a  set  of  five  diffusions 
carried  out  at  temperatures  from  650-750  *C.  The  plots  have 
a  normalized  depth  scale  to  show  the  change  in  profile  shape. 
To  do  this,  we  defined  jc,  as  the  point  on  each  curve  where  the 
Zn  concentration  equals  1019  cm-3,  then  divided  the  depth 
scale  for  each  curve  by  this  xr  The  relative  size  of  the  tail 
region  grows  larger  at  higher  diffusion  temperatures.  The 
650  °C  diffusion  is  the  most  desirable  from  a  device  process¬ 
ing  standpoint,  since  it  has  the  highest  surface  concentration 
and  a  small  tail  region.  Hole  concentration  for  the  650  ”C 
diffusion  was  1.3 X  10:o cm-3. 

Figure  3  is  a  set  of  SIMS  profiles  for  diffusions  earned 
out  at  700  °C  for  various  times  from  1-10  min.  Figure  4 
shows  the  progression  of  three  points  on  each  profile  plotted 
versus  the  square  root  of  time.  This  figure  suggests  that  the 
entire  Zn-diffusion  profile  scales  with  the  square  root  of  dif¬ 
fusion  time,  so  that  the  general  profile  shape  does  not 
change.  This  is  further  suggested  by  Fig.  5,  in  which  the  four 
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FIO  2.  SIMS  profiles  for  a  set  of  live  diffusions  earned  out  on  era  tempera¬ 
ture  ranee  from  b50-750°C  -  1-0  Torr.  /»ol/ll  -  :  0  Torr  r-.e 

piois  ha\e  a  normalized  depth  scale  to  show  the  chance  in  profile  shape  Th; 
l.x/.x  i  =  1  point  on  this  scale  is  the  place  where  the  Zn  concentration  s 
10'  ‘cm  the  actual  depth  to  which  tnis  corresponds  is  shown  tn  par.r.’he- 
-  is 
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FIG.  3.  SIMS  profiles  for  a  set  of  four  diffusions  carried  out  at  700  ‘C  for 
times  from  1-10  min. 

profiles  of  Fig.  3  are  replotted  on  a  normalized  depth  scale. 
The  error  bars  in  Fig.  5  show  the  maximum  deviation  of  any 
profile  from  the  average  of  the  four.  Within  our  experimen¬ 
tal  error,  the  shape  of  the  diffusion  profile  does  not  change 
with  time.  Dohsen  et  al.2  show  a  similar  result  to  that  in  Fig. 
4.  Figures  3-5  also  indicate  that  relatively  little  movement  of 
the  Zn  profile  occurs  during  the  1-min  purge  after  the  DEZn 
flow  is  switched  off.  If  significant  movement  did  occur,  one 
would  expect  the  relative  change  to  be  much  greater  in  a  1- 
min  diffusion  than  a  10-min  one. 

For  the  purpose  of  comparing  the  various  diffusion  pro¬ 
files,  it  is  convenient  to  define  a  “kink  level"  where  the  two 
different  portions  (the  surface  region  and  the  tail)  of  each 
profile  meet.  Each  profile  has  a  small  section  in  which  it  is 
concave  up;  the  kink  level  is  defined  as  the  point  where  the 
second  derivative  is  at  a  maximum.  The  “x”  in  Fig.  3  is  an 
example  of  the  kink  level.  Figures  6  and  7  show  the  kink  level 
of  several  profiles  plotted  against  diffusion  time  and  tem¬ 
perature.  Kink  level  rises  with  temperature  and  is  approxi- 
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FIG.  5.  The  four  SIMS  profiles  from  Fig.  3  are  shown  plotted  on  a  normal¬ 
ized  depth  scale.  The  depth  at  which  the  Zn  concentration  of  each  profile 
falls  to  10|7  cm-3  is  set  equal  to  1.0.  The  error  bars  show  the  maximum 
deviation  of  any  profile  from  ihe  average  of  the  four  The  change  in  profile 
shape  with  time  appears  to  be  very  small. 


mately  independent  of  time.  Figure  7  also  shows  the  vari¬ 
ation  of  surface  concentration  with  temperature. 

We  have  studied  the  effect  of  DEZn  and  TMAs  partial 
pressures  on  the  diffusion  profiles.  For  a  range  ofpDEZ„  from 
0.2-2.0  Torr,  Zn  surface  concentration  increases  linearly 
with  /tDEZ„-  If  we  ignore  the  tail  region  of  the  profiles  and 
apply  the  interstitial-substitutional  (I-S)  model  to  surface 
region  (C>  10,Q  cm-3),  we  can  extract  a  surface  diffusivity 
according  to  the  equations  of  Weisberg  and  Blanc4: 
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FIG  4  The  junction  depth  of  various  Zn  concentrations  (  V  =  I01’;  FIG.  0  Kink  levelasa  function  ofdiffusion  time.  Kink  level  isdelined  asthe 

.V  =  10'".  and  A  =  lO'"  cm  ')  for  the  SIMS  profiles  shown  in  Fig  3.  plot-  point  where  the  two  different  pomons  ( the  surface  region  and  ihe  tail)  of 

led  vs  the  square  tool  of  diffusion  lime.  each  profile  meet.  An  example  is  the  ".t"  in  Fig  3 
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FIG.  7.  Kink  levei  and  surface  concentration  as  a  function  of  inverse  diffu¬ 
sion  temperature. 

Au,  =  <x,/1.092)2(l/r), 

D  =  Dm{C/Cm)\ 

Figure  8  plots  this  extracted  surface  diffusivity  against  the 
DEZn  partial  pressure.  The  slope  of  the  straight  line  is  2, 
indicating  a  good  fit  to  the  interstitial-substitutional  model 
in  the  surface  region.  TM  As  partial  pressure  has  only  a  small 
effect  on  the  profiles  over  a  range  from  0.2-3.0  Torr.  At 
650  °C,  Aur  decreases  as  approximately  the  fourth  root  of 
p rMAl  as  the  I-S  model  predicts. 
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MG  8  Surface  diffusivity  (Au»  )  plotted  against  DEZn  partial  picssure. 
Surface  diffusivity  is  defined  by  the  equation  D  lf  —  f  v  / 1  002  )■  ( I//).  To 
D  ,.  <  is  taken  as  I  he  point  where  C-- ( Zn  |  .=  I0,gcm 


Figure  9  shows  the  variation  of  kink  level  with  Zn  sur¬ 
face  concentration  (Csur).  C„ur  varied  by  adjusting  the 
DEZn  partial  pressure.  The  kink  level  varies  as  approxi¬ 
mately  the  0.67  power  of  Csur,  although  the  power  depen¬ 
dence  could  actually  be  from  0.5  to  1.0  if  experimental  error 
is  considered. 

IV.  MODEL 

Double  profile  or  kinked  Zn  diffusions  have  been  ob¬ 
served  by  Tuck  and  Kadhim.5  Their  data  also  indicate  a 
change  in  profile  shape  with  diffusion  time,  which  they  at¬ 
tribute  to  the  gallium  vacancy  concentration  falling  below  its 
thermal  equilibrium  value.  Similar  data  and  a  more  general 
model  are  discussed  in  a  paper  by  Ting  and  Pearson.8  Gosele 
and  Morehead1'  propose  a  slightly  different  mechanism  to 
explain  the  data  of  Tuck.  They  suggest  that  the  dominant 
defect  in  GaAs  may  be  interstitial  gallium  (/G, )  and  that  the 
kink  occurs  because  /G,  is  present  in  excess  of  its  thermal 
equilibrium  value  both  near  the  diffusion  front  and  in  the  tail 
region.  These  models  are  all  quite  plausible,  and  the  non- 
equilibrium  they  describe  probably  does  occur  under  certain 
conditions.  However,  the  concentration  dependence  of  the 
diffusion  coefficients  they  predict  does  not  exactly  match 
our  data  (or  Tuck's).  We  propose  that  the  double  profile  is 
caused  by  the  existence  of  two  different  species  of  diffusing 
Zn  (in  addition  to  substitutional  Zn)  which  dominate  in 
different  locations  according  to  the  position  of  the  Fermi 
level. 

We  begin  our  analysis  by  determining  the  effective  Zn 
diffusion  coefficient  as  a  function  of  concentration  from  the 
SIMS  profiles.  Tuck  and  Kadhim5  have  pointed  out  that  the 
standard  Boltzmann-Matano  analysis10  is  not  valid  for  a 
profile  shape  which  changes  with  time  Since  our  data  do  not 
indicate  the  time  dependence  that  Tuck  observed,  we  have 
applied  the  Boltzmann-Matano  method,  but  it  is  probably 
best  to  regard  the  result  as  an  approximation.  Our  data  do 
not  prove  that  there  is  no  time  dependence,  only  that  it  is 
very  small  under  our  experimental  conditions. 

Figure  10  gives  the  result  of  our  Boltzmann-Matano 
analysis.  In  the  surface  region  we  have  D  -  ~  C  ~  as  predicted 
by  the  I-S  model  There  is  a  transition  region  between  about 
2  ■  inland  2  ■  1 0 1  '*  cm  'in  which  A..  -  C'  -.  The  tail  re- 
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FIG.  10.  Effective  Zn-diffusion  coefficient  as  a  function  of  concentration  for 
the  SIMS  profile  in  Fig.  1,  obtained  by  Boltzmann-Matano  analysis. 

gion  is  noisy  and  it  is  difficult  to  determine  whether  we  have 
Dts  “Cor  Drfr  cC1,  although  for  the  particular  data  in  Fig. 
10,  oc  C  appears  to  give  a  better  fit.  The  concentration 
dependence  can  also  be  determined  by  fitting  calculated  dif¬ 
fusion  profiles  to  the  SIMS  data.  The  profiles  which  result 
from  assuming  either  oc  C 2  or  Air«C  can  be  calculated 

by  the  method  of  Weisberg  and  Blanc.4  Figure  1 1  shows  the 
SIMS  data  of  Fig.  1  fit  with  a  Detl  cr  C1  profile  in  the  surface 
region  and  a  Dt„  <x  C  profile  in  the  tail.  For  all  of  our  SIMS 
data,  assuming  a  Dc„  oc  C  profile  gives  a  much  better  fit  to  the 
tail  region  than  De„  <x  C2.  DtK  oc  C also  appears  to  give  a  bet¬ 
ter  fit  to  the  tail  region  of  Tuck’s  data,  shown  in  Fig.  4  of 
Tuck  and  Kadhim.5  Neither  Tuck’s  nor  Gosele’s  model  for 
Zn  diffusion  predicts  this  change  in  concentration  depen¬ 
dence. 

In  the  I-S  model,  substitutional  Zn  occupies  the  Ga  sub- 


FIG  1 1  Calculated  diffusion  profiles  have  been  tit  to  the  SIMS  data  of  Fir 
1  Assuming  a  u  C:  gives  a  good  til  to  the  surface  region:  D.,  a  C 

better  in  the  tail 


lattice  as  a  shallow  acceptor  (Zn,  ),  while  interstitial  Zn  is 
considered  to  be  a  singly  or  doubly  ionized  donor  (Zn,'  or 
Zn,'  *  ).  The  solubility  of  Zn,“  is  much  greater  than  that  of 
interstitial  Zn.  while  interstitial  Zn  has  a  diffusion  coefficient 
several  orders  of  magnitude  greater.  Interstitial  Zn  enters 
the  lattice  at  the  surface,  moves  rapidly  to  the  diffusion  front, 
and  becomes  substitutional  according  to  either 

Zn,  -r  I  g»  ^  Zn,  -i -  t.h 
or 

Zn,4  4  +  Vg,^  Zn  “  +  Ih  4, 

where  h  4  are  holes.  Both  of  these  reactions  yield  Dt„  <x  C:. 
Precisely  which  reaction  occurs  is  not  agreed  upon.  Seem¬ 
ingly  contradictory  evidence  (from  isoconcentration  stud¬ 
ies)  exists  for  both  Zn,4  and  Zn,4  4  .“-12 

Figure  1 1  tends  to  suggest  that  two  different  diffusing 
species  or  mechanisms  are  involved  in  the  diffusion.  We  pro¬ 
pose  that  Zn  can  exist  as  either  a  singly  ionized  interstitial 
(Zn,4  )  or  a  doubly  ionized  interstitial  (Zn,4  4  ),  in  addition 
to  being  a  substitutional  acceptor.  Our  model  assumes  that 
Zn,  is  a  deep  donor  level,  within  a  few  tenths  of  an  eV 
from  the  valence-band  edge.  Zn,*  would  also  be  a  deep  do¬ 
nor,  but  near  the  center  of  the  band  gap.  We  make  this  choice 
of  donor  level  based  on  a  study  of  Zn  diffusion  in  heavily  n- 
type  GaAs  by  Wang.1-’  Wang  observed  only  substitutional 
Zn  diffusion  and  concluded  that  Zn,'  does  not  exist  in  n- 
tvpe  material.  This  places  the  Zn,'  donor  level  at  kast  a  few 
tenths  of  an  eV  from  the  conduction-band  edge,  and  prob¬ 
ably  deeper.  The  energy  difference  between  Zn, 4  andZn  *  ' 
can  be  estimated  by  the  method  of  Harrison,14  taking  the 
difference  between  the  first  and  second  vacuum  ionization 
potentials  for  Zn  and  dividing  by  the  dielectric  constant  for 
GaAs: 

±E^(E2-E,)/eatAl  =0.65  eV. 

These  donor  levels  are  illustrated  in  the  band  diagram  of  Fig. 
12. 

The  Zn*  4  would  be  formed  at  the  surface  and  domi¬ 
nate  the  diffusion  process  in  the  heavily  doped  surface  re¬ 
gion.  In  less  heavily  p- type  regions  near  the  diffusion  front, 
the  Fermi  level  would  move  away  from  the  valence-band 
edge.  The  more  mobile  Zn,4  would  then  be  favored,  giving 
rise  to  the  diffusion  tail  and  yielding  a  D^„  a  C  dependence. 
This  is  summarized  in  the  following  reactions.  We  assume 
that  interstitial  Ga  is  the  dominant  defect,  after  Gosele.3 
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FIG  12.  GaAs  band  diagram,  showing  ihe  position  of  the  two  donor  levels 
assumed  in  the  model  The  E  '  '  level  for  Zn  is  choosen  to  be  about  0  2 
eV  from  the  v alenee-band  edge  as  the  best  fit  to  our  data.  The  E  *  level  for 
Zn  *  ts  estimated  to  lie  about  0  ft?  eV  higher  in  energy,  based  on  relative 
iv-ntzation  potential 
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These  reactions  would  also  work  if  we  assumed  V  G,  to  be 
the  relevant  defect: 

Zn/  +  —  Zn/  +  /q«  +  2h*,  (la) 

Zn/  +  —  Zn/  +h*,  (2a) 

Zn/  —  Zn,"  +  Iq»  +  h  *■  (3a) 

Note  that  Zn/  +  must  have  a  smaller  diffusion  coefficient 
than  Zn,'*'  in  order  for  the  model  to  work.  This  is  physically 
plausible  if  the  main  barrier  to  interstitial  motion  is  electro¬ 
static.  Although  Zn/  +  has  a  smaller  radius  than  Zn/,  they 
are  probably  both  small  compared  to  gaps  in  the  GaAs  lat¬ 
tice.  There  is  the  possibility  that  doubly  ionized  Zn  forms 
some  sort  of  charged  complex  with  Ga  or  As  atoms  on  the 
GaAs  lattice,  or  with  a  point  defect.  The  notim  of  a 
(Zn,  —  Ga)  +  +  or  (Zn,  —  As)  +  +  complex  is  one  way  to 
conceptualize  an  electrostatic  barrier  to  interstitial  motion. 
Our  model  does  not  require  the  existence  of  any  such  com¬ 
plex,  but  only  that  the  doubly  ionized  Zn  diffuse  more  slowly 
than  Zn/ . 

The  concentration  of  Ia ,  is  assumed  to  be  everywhere 
equal  to  its  equilibrium  value.  The  following  equations  then 
hold  if  reactions  ( la)-(3a)  are  in  equilibrium: 

cu=k,cl  ;ib) 

Cu  =  k2C,C„  (2b) 

C,=kyC],  (3b) 

where  C„  =  [Zn/  *  ],  C,  =  [Zn/  ],  and  C,  =  [Zn/  ]. 
The  parameters  kx  through  k,  are  equilibrium  constants. 
The  observation  that  Zri  diffuses  rapidly  only  when  it  is  be¬ 
ing  incorporated  at  the  surface  implies  that  reactions  ( la)- 
(3a)  go  predominantly  in  one  direction,  as  indicated  by  the 
large  arrows.  This  is  intuitively  plausible  because,  in  each 
case,  going  in  the  opposite  direction  would  appear  to  require 
a  higher  activation  energy.  For  instance,  Zn  interstitials 
might  be  readily  formed  at  the  surface  because  of  unsatisfied 
bonds  and  a  large  excess  of  Zn-  Once  in  the  lattice,  formation 
of  Zn/  +  or  Zn/  from  Zn/  would  require  breaking  bonds 
to  the  neighboring  Ga  and  As  atoms  and  leave  a  vacancy. 

Under  these  conditions,  our  computer  simulation  is  nec¬ 
essarily  an  approximation  to  the  behavior  of  reactions  ( la)- 
(3a).  Our  approximation  is  that  the  deviations  of  C„  or  C, 
above  the  equilibrium  concentrations  predicted  by  ( lb)— 
(3b)  are  corrected  immediately  via  reactions  (la)-(3a). 
Deviations  below  the  equilibrium  value  are  allowed  to  per¬ 
sist  (although  they  are  corrected  by  diffusion).  We  solve  the 
following  equations  for  simultaneous  diffusion  and  chemical 
reaction: 


dC„ 

dt 

<k 

dt 

3C,  dS i 
dt  dt 


=  Z>„ 


=  D, 


dlC 

w  ri 

dx: 

d2C, 


di r, 

’  dt 

ds? 

dt  ' 


dS , 

'  dt’ 
dS, 
dt' 


(4) 


(5) 


(6) 


dx~ 
dS, 
dt' 

D„  and  D,  are  the  diffusion  coefficients  for  doubly  and  singly 
ionized  Zn  interstitials.  (dSt/dt)  is  the  rate  at  which  Zn/  * 
become  Zn/  according  to  reaction  ( la).  dSWdt  is  the  rate  at 


FIG.  13.  Results  of  a  computer  simulation  for  a  5-min  diffusion  at  725  *C, 
showing  the  concentrations  of  Zn~ ,  Zn,"  * ,  and  Zn,"  Note  that  the  total 
concentration  of  mobile  Zn,  which  is  assumed  in  the  model  to  be  much  less 
than  [  Zn,"  ],  is  arbitrarily  taken  to  be  10“ 3  of  the  total  Zn  present  at  the 
surface. 
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FIG.  14  SIMS  data  and  a  fitted  computer  simulation  for  a  5-min  diffusion 
at  725  -C. 
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TABLE  I.  The  model  parameiers  used  10  fit  our  experimental  data  of  Fig  2. 


Model 

parameter 

Units 

Temperature  (’C) 

650 

675 

700 

725 

750 

C„, 

cm  ' 

2.4  X  10:o 

1.8x10-’" 

I.3X10-" 

1.3  X  10;" 

6 

X 

«N 

D„ 

cm;/s 

1.2x  10’” 

0.8X10-” 

0.4  X  10“ 12 

1.2  X  10-'-’ 

0.7  X  10- i: 

D 

cm2  /r» 

-■'Ox’"-'- 

3.0 X  I0-|! 

1.5X10” 

5.0x10-'-’ 

6  Ox  10  ~ 

Ik-)-' 

cm  1 

1  IX  10'“ 

1.9x10'“ 

3. OX  10'* 

4  5x  10" 

7. Ox  10" 

)  1 

no  unit 

20 

17 

15 

12 

9 

which  Zn*  +  become  Zn,+  according  to  reaction  (2a).  Fin¬ 
ally,  dSffdt  is  the  rate  at  which  Zn,*  become  Zn  "  via  reac¬ 
tion  (3a). 

Equations  (4)-(6)  are  solved  by  a  finite  difference  algo¬ 
rithm.  selecting  small  values  for  At  and  Ax.  For  each  small 
time  step  At,  the  mobile  species  are  allowed  to  diffuse  ac¬ 
cording  to 


AC„  (x„ )  =  At  (d, 
AC,  (x„ )  =  At  (A 


dx:  )' 

OiCJxm)\ 

dx'-  J 


Next,  reaction  (2a)  is  assumed  to  bring  C„  (x„ )  and  C,  (x„ ) 
into  equilibrium.  AS2(x„ )  is  obtained  by  calculating  the  ex¬ 
pected  equilibrium  value  for  C'q  (x„ )  from  Eq.  (2b),  then 
taking  the  difference  A5:(x, )  =  C„  (x„ )  —  C'q  (x„ ).  A S2 
is  constrained  to  be  >0,  since  we  are  allowing  reactions 
( la)— ( 3a)  to  proceed  only  in  one  direction.  C„(x„)  and 
C,(. x„)  are  adjusted  using  AS-.  Next,  reactions  (la)  and 
(3a)  are  assumed  to  bring  C„(x, )  and  C,  (x„)  each  into 
equilibrium  with  C,  (x, ).  AS,(x„)  and  A S,(x„)  are  ob¬ 
tained  as  AS:(xn)  was  before.  Finally,  C„  (x„ ),  C,  (x„  ), 
and  C,  (x„  )  are  updated  and  t  is  incremented  by  Ar. 

Toobtain  initial  conditions,  reactions  ( la)  and  (2a)  are 
assumed  to  be  in  equilibrium  exactly  at  the  GaAs  surface. 
The  Zn  surface  concentration  C.  (x  =0)sC,u[  is  known. 
Choosing  A,  fixes  C„  (x  =  0)  and  C,  (x  =  0)  within  an  arbi¬ 
trary  constant.  The  other  variable  parameters  are  A,,  D  , , 
and  D,  (for  a  total  of  four  at  each  temperature).  Figure  13 
shows  the  result  of  a  computer  simulation  for  a  5-min  diffu¬ 
sion  at  725  °C,  including  the  concentrations  of 
Zn,' .  Zn,*  * ,  and  Zn/  .  Figure  14  shows  SIMS  data  for  a  5- 
min  diffusion  at  725  *C  along  with  a  computer  calculated 
profile. 

This  model  has  been  used  to  fit  our  experimental  data 
over  a  temperature  range  from  650  to  750  °C,  with  results 
similar  to  Fig.  14.  For  any  given  temperature,  this  model 
predicts  a  constant  profile  shape  which  progresses  into  the 
GaAs  as  t l/:,  in  accordance  with  our  data.  Table  I  lists  the 
model  parameters  used  to  fit  each  temperature.  D.,  and  D ,  do 
not  show  any  definite  trend  with  temperature,  reflecting  the 
data  of  Fig.  ..  This  observation  h3S  been  made  before  with 
Zn  diffusions  under  other  conditions.  "-M  and  is  often  asso¬ 
ciated  with  an  interstitial  diffusion  mechanism.  Note  that  D 


and  D,  are  effective  diffusivities,  meaning  that  they  depend 
on  both  the  interstitial  diffusion  rate  and  the  interstitial  con¬ 
centration.  The  erratic  behavior  of  D„  and  D,  with  tempera¬ 
ture  is  probably  due  to  a  decrease  in  interstitial  concentra¬ 
tion  with  temperature  (since  Zn  surface  concentration  falls 
as  temperature  rises)  and  an  increase  in  interstitial  diffusi- 
vity.  The  value  of  parameter  k2  is  consistent  with  assuming 
that  Zn,*  *  is  a  deep  donor  roughly  0  7  eV  from  the  valence 
edge,  the  exact  donor  level  varying  with  temperature.  Since 
our  model  has  four  parameters  which  affect  profile  shape, 
fitting  the  data  with  these  four  parameters  is  not  by  itself  an 
adequate  test  of  the  model.  A  good  test  would  be  an  indepen¬ 
dent  determination  of  the  effective  diffusion  coefficients  D„ 
and  D„  especially  their  variation  with  temperature.  More 
experiments  involving  the  variation  of  the  proh.e  shape  with 
time  would  also  be  valuable. 

V.  SUMMARY  AND  CONCLUSIONS 

We  have  studied  the  diffusion  of  Zn  into  GaAs  using 
diethylzinc  and  trimethylarsenic  sources.  This  method 
yields  good  control  over  Zn  surface  concentration  and  junc¬ 
tion  depth.  The  variation  of  the  diffusion  profile  with  time, 
temperature,  and  alkyl  partial  pressures  has  been  examined. 

These  diffusions  have  a  complex  double  profile  which  is 
not  predicted  by  the  interstitial-substitutional  model  for  Zn 
diffusion.  We  suggest  that  this  double  profile  may  be  caused 
by  the  existence  of  two  mobile  Zn  species  which  dominate  in 
different  regions  according  to  the  position  of  the  Fermi  level. 
We  have  implemented  a  computer  simulation  which  gives  an 
approximate  solution  based  on  this  two-species  model.  The 
simulation  gives  a  good  fit  to  our  experimental  data. 

Zn  diffusion  in  GaAs  is  a  complicated  process.  Present 
experimental  data  do  not  completely  validate  our  model  nor 
the  existing  alternatives.  The  lack  of  time  dependence  (of  the 
profile  shape)  we  observe  may  simply  be  a  special  case  of 
Tuck's5  or  Gosele's9  model  for  a  short  diffusion  times.  Our 
model  shows,  however,  that  a  double  diffusion  profile  does 
not  necessarily  imply  vacancy  (or  interstital)  nonequilibri- 
um. 
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